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Iwc methods for evaluating the accuracy of hydrographic 
positioning data are presented. One method consists of 
classifying each position in a survey based on the radius of 
the 90 percent confidence circle. The second method 
involves classification of positions based on the parameters 
of the 90 percent confidence ellipse. Both methods are 
based on geometric and statistical relationships tetween 
intersecting lines of position. 

Bange-range, azimuth-azimuth, and range-azimuth posi- 
tioning data are classified using both criteria. For 
noncritical positions, the confidence circle method is found 
to be preferable due to its ease of interpretation. For 
positions of significant features, such as underwater 
hazards, the confidence ellipse provides a more useful 
representation of the shape and orientation of the true 
error distribution. 

The concept of presurvey positioning design is alsc 
presented. With the aid of computer graphic displays, the 
hydrographer can predict the accuracy of offshore posi- 
tioning data prior to data acquisition. By analyzing accu- 
racy lobes generated about shore stations, a survey can be 
designed to meet given specifications. 
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I. INTEODDC^ON 



a. EACKGRCOND 

A hydrographic record can he viewed as the resultant of 
two independent measurements made at a discrete point over a 
body of water. These measurements involve the determination 
of a vessel’s position at a given time as well as the depth 
of water at that position. Of interest to the hydrographer 
and to the user of hydrographic data is the accuracy of the 
position determinations. Fundamental to the determination 
of positional accuracy is the identification of the sources 
of errors in position measurements and the ultimate treat- 
ment of these errors. 

A hydrographic position can be determined by a number of 
methods all involving geometric relationships between known 
points and the vessel’s unknown location. The known points 
may be fixed stations on shore, whose coordinates have been 
determined by geodetic survey methods, or they may be 
rapidly moving satellites whose coordinates in time and 
space can he defined very precisely. A hydrographic posi- 
tion is established by the intersection of two or more lines 
of position (LOP’s) which are generated by the geometric 
relationships between the fixed points and the vessel’s 
unknown location. The resultant accuracy of the vessel’s 
position is therefore, in part, a function of the errors 
associated with the irtersecting LOP’s. 

Several measures cf accuracy can be used to evaluate the 
guality of a hydrographic position. Predictability, or 
absolute accuracy, is the measure of accuracy with which the 
positioning system can define the location of the same point 
in terms of geographic coordinates. Repeatability, or 
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relative accuracy, is a measure with which a positioring 
system permits a user to return to a specific point on the 
earth’s surface in terms of the LOP’s generated fcy the 
system [Bef. 1, p. 1 ]. With the elimination of all system- 
atic cr tias errors, the terms repeatability and predicta- 
bility become identical. Hydrographic surveyors usually 
work toward this condition, although it is not always 
achievable. 

Heinzen [Ref. 2] and Burt [Bef. 3] have presented 
several techniques for quantifying the repeatable accuracy 
for offshore positions. These techniques have roots in tbe 
statistical treatment of random error. Although the methods 
have teen well documented, no single criterion to classify 
the accuracy of a hydrographic position has been agreed upon 
by the international hydrographic community. 

Preceding the development of automation in hydrographic 
data acquisition and processing, the task of calculatirg an 
accuracy figure to attach to each position in a hydrographic 
survey was unthinkable. To ensure overall accuracy in a 
survey, certain generalizations were developed to act as 
guidelines. For example, the U.S. Coast and Geodetic Survey 
Hydr ograrhic M an u al [Bef. 4, p. 217] states the following 
concerning the strength of a three-point fix: 

The fix is strong when the sum of the two angles is 
equal to or greater than 18Qo and neither angle is less 
than 3C0. The nearer the angles equal each other the 
stronger will be tbe fix. 

Generalizations of this type provided useful qualitative 
guidance for assuring a degree of positional accuracy and 
many are still in existence today. 

With the aid of computers, the hydrographer now has the 
capacity to evaluate the accuracy of positioning data fcr an 
entire survey. An accuracy figure can be computed for each 
position in a survey and stored in a data base along with 
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ether survey inf ormation. This figure may provide useful 
inf or HI at ion for users of the data, as well as a yardstick 
for the hydrographer to evaluate the quality of the work. 
Furt her mere, a presurvey accuracy analysis enables a survey 
to be designed to meet desired specifications. 



E. ACCOEACY STANDABES FOB HYDBOGBAPHIC POSITIONING 



In 1S82, the International Hydrographic Organization 
(IHO) putlished new recommendations for error standards 
concerning the accuracy of hydrographic positions. These 
standards [Bef. 5] are; 



The position of soundings, dangers and all other signif- 
icant features should be determined with an accuracy 
such that any probable error, measured relative to shore 
contrcl, shall seldom exceed twice the minimum plottatle 
error at the scale cf the survey (normally 1.0 mm on 
paper) . It is most desireable that whenever positions 
are determined by the intersection of lines of position, 
three such lines be used. The angle between any pair 
should not be less than 30°. 



Most statisticians define the term "probable errer" as 
that error occurring at the 50 percent probability level. 
However, the author cf the IHO standards. Commodore A.H. 
Cooper BAN (Ret.) has stated that the term "probable error" 
was intended to have no statistical significance. Munson 
interpreted the words "shall seldom exceed" to mean 10 
percent of the time £Eef. 6]. Using this interpretation, 
the first sentence of the specification might be written; 



The position of soundings, dangers and ail other signif- 
icant features should be determined with an accuracy 
such that any error in position measured relative to 
shore control will fall within a circle with radius of 
the minimum plottable error at the scale of the survey 
(normally 1.() mm. on paper), with 90 percent confidence. 



The specification in this form could be evaluated quantita- 
tively. The criterion for defining accuracy in terms of a 
fixed probability is common in the field of surveying. For 
example, the standards of accuracy developed for geodetic 



contrcl surveys have their origin in probability theory. 
Procedures for obtaining first-order geodetic positions 
require sixteen repeated theodolite observations of each 
direction. Lower order positions require fewer numbers of 
observations. Given the precision of one observation of 
each direction, it can be demonstrated that increasing the 
number of observations coincides with increasing the prob- 
ability of the direction falling within specified limits. 

Eegarding accuracy determinations, there are several 
problems unique to hydrographic surveying. Whereas stan- 
dards for other types of surveys rely on multiple observa- 
tions of the same quantity, the accuracy of a hydrographic 
position must be evaluated in terms of a single observation 
(which may be the intersection of two or more LOP’s). 

Diverse methods for obtaining a hydrographic position exist 
and these methods must all be evaluated using the same 
criterion. Also, there is a broad spectrum of equipment 
used in hydrographic positioning and in many cases the 
precision of this equipment is not well defined. 

C. OBJECTIVES 

A need exists to give quantitative meaning to the accu- 
racy specifications set forth by the IHO. One of the objec- 
tives of this thesis is to demonstrate that defining the 
specifications in terms of the fixed 90 percent confidence 
level is a valid interpretation. By defining what the spec- 
ifications imply, procedures can be developed to meet the 
stand ards. 

A second objective of this thesis is to apply the theory 
of errors, associated with hydrographic positioning, to a 
data set. This analysis involves classifying positioning 
data acquired in a survey based on the radii of circles of 
equivalent probability. It will be demonstrated that this 
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coethod of classification is a useful index for quantifying 
the accuracy of positions. The computed radii of the 90 
percent confidence circles can serve as an accuracy figure 
that can be attached to each position in a survey and stored 
in a data base. 

The third objective of this thesis is to demonstrate 
that a presurvey analysis can be used in designing posi- 
tional accuracy to meet specifications. The existing 
general guidelines for planning can be tetter defined. For 
example, in planning a survey hydrographers usually lay out 
circles Vihich delimit the 3Qo and 150° boundaries that 
define the minimum and maximum allowable intersection angles 
between two LOP’s. As a means to meet accuracy require- 
ments, it can be shown that these limits should vary based 
on the scale of the survey and the precision of the posi- 
tioning equipment. 
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II. MI5M Ql the p bobl em 



The development cf an accuracy figure for offshore posi- 
tions is inherently tied to the geometry of the positioning 
method and the errors which are associated with the posi- 
tioning eguipment that is used. This chapter will discijss 
the geometric and statistical elements involved in deter- 
mining an offshore position and presents several methods for 
quantifying repeatable accuracy. 

A. HYDECGEAPHIC POSITIONING GEOMETRIES 

An offshore fix can be determined by the intersection cf 
two or more LOP’s. These LOP's may be generated by elec- 
tronic or visual means. Working toward the development of 
an accuracy index, it will be necessary to compute the angle 
of intersection of the LOP's associated with different posi- 
tioning geometries. The following sections discuss the 
geometry cf conventional offshore positioning methods and 
ways to compute the argles of intersection. This thesis 
will not address the geometry involved in a three-point 
sextant fix, 

1 , Ra n ge-Range 

Establishing an offshore fix by range-range geometry 
involves measuring distances electronically from fixed posi- 
tions on shore to the vessel's unknown location. Ranges can 
he determined by measuring the elapsed time between trans- 
mission and receipt of a radio pulse or by comparing the 
phase of the transmitted wave with the phase of the received 
wave [Ref. 2]. In each case, transmitters are set on 
stations on shore whose coordinates are determined by 
precise land survey methods. 



An electronic positioning system may be active or 
passive. In an active system, a transmitter from the survey 
launch keys the transmission of ranges from the shore 
station. In turn, the signals generated from the shore 
stations (slaves) are then received by the launch. An 
active system is limited to a finite number of users, 
usually not more than about four. The number of users of a 
passive system is unlimited as the survey launch requires 
only a receiver which is constantly listening for signals 
which are being transmitted from shore. 

Short-range, or line-of-sight, positioning systems 
are used for nearshore hydrographic surveys. These systems 
operate in the microwave region of the electromagnetic spec- 
trum (3 to 10 GHz). A distance is determined by observing 
the time needed for a pulse to travel from a master tran- 
sponder located aboard the survey vessel to a remote tran- 
sponder on shore and back to the master transponder. 

Knowing the average velocity of the electromagnetic pulse, 
the distance D is then 




( 2 . 1 ) 



where c is the group velocity of the wave packet and t is 
the two-way travel time. Short-range systems which are in 
wide use today are Bacal Decca’s "Trisponder" and Motorola’s 
"Mini-Hanger.” These systems have direct range readout and 
are readily interfaced into a navigational computer and a 
data acquisition system. Both systems are active and user 
limited . 

Medium-range positioning systems operate in the 1- 
to 5-KHz frequency range of the electromagnetic spectrum. A 
distance is determined by measuring the phase relationship 
between transmitted and received waves. These systems are 
usually referred to as continuous wave systems and the 
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problem cf lane ambigoity must be addressed. Ranges are 
expressed in full and partial lane counts where a lane width 
w is 



w 




( 2 . 2 ) 



where X is the wavelength of the transmitting frequency, f, 
and given by 



X = (2.3) 

Medium-range systems ccmmonly in use today are Cubic 
Western’s "ARGO," Hasting Raydist’s ••Raydist," and Odom 
Offshore’s ”Hy drotracX. ” 

Ihe angle of intersection associated with a range- 
range position is computed frcm a simple trigonometric rela- 
tionship. The vessel’s position P (Fig. 2.1) is determined 
by the intersection of the ranges from the left and right 
shore stations, R1 and R2 respectively. 3 is the base line 
distance computed between the two known shore stations. 

Since the range circles from the shore stations intersect at 
two points, it is necessary for the plotter to recognize 
which side of the base line the vessel is on in order to 
eliminate the ambiguity. The angle of intersection cf the 
two LOP’S (B) is given by the law of cosines 

6 » 180° - Arc cos ( ^ (2.4) 

In qualitative terms, the fix is strongest when B approaches 
900. Host hydrographic specif ications limit the angle cf 
intersection from a minimum of 30 ° to a maximum of 150®. 
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Figure 2. 1 Gecnetry of a Eange-Hange Position 

2 • H yperb o lie- H vee r boli c 

Hydrographic positioning by hyperbolic-hyperbolic 
geometry utilizes the intersection of two hyperbolas each 
generated about a pair of shore control stations. A hyper- 
bola is the locus of points in which the difference of 
distance from two fixed points is always constant. A three- 
station hyperbolic net is the most commonly used hyperbolic 
mode for offshore survey (Fig. 2.2). One family of hyper- 
bolas (Red) are generated about a master station, M, and a 
slave, 3; while a second family of hyperbolas (Green) are 
generated with respect to the master and a second slave, G. 
For the first family cf hyperbolas, the control points M and 
R act as the foci, while points M and G act as the foci for 
the second family. 
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Hyperbolic location methods can be divided into tvo 
groups based on the electronic principles used to define the 
distance differences £Eef, 1 , p. 87]. Loran is an example 
of a pulse system in which the differences in times cf 
arrival of pulses transmitted by the master-slave combina- 
tions are translated into distance differences. The resul- 
tant position has no lane ambiguity and is easily resolved. 
The second method of hyperbolic positioning involves meas- 
uring a phase difference from two master-slave combinations 
at the vessel's position. The phase difference translates 
into a fractional lane count which in itself provides an 
ambiguous position. This ambiguity is resolved by using a 
whole-lane counter which is initialized at a known geograph- 
ical point. In hyperbolic positioning, the ship is in a 
passive mode and the system can be used by many vessels. 

The angle of intersection between the two hypertolas 
can be computed by first defining the following quantities: 

S is the length of red base line, 

S is the length of green base line, 

R is the distance between master and vessel's pcsiticn ?, 
is the distance from red slave to point P, 

R is the distance from green slave to point P, 

o 

is the angle between lines PM and PR, and 

a is the angle between lines PM and PG. 
g 

The spacing between lanes increases with distance 
from the master-slave pair. The lane widths along the base 
line are 



w' 

r 



and 



w' X 

3 V 



(2.5) 



Then the lane widths at any point P are 

/ 1 \ and ^ ^ 1 ) 

IT ^ sin (a /2) ^ ^ 2 ^ sin (a /2) (2.6) 



w 
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■where the term 1/sin ( a/2) is called the lane expansion 
factor. The angle of intersection ,8, between the two 
hyperbolas is then given by 



6 = 



o.. 

r q 
2 



(2.7) 




Figure 2.2 Geometry of a Hyperbolic-Hyperbolic Position 



5 . Banq e- Azimu th 

This positioning geometry is used for nearshore, 
line-cf- sight surveys. One LOP is generated by an elec- 
tronic range originating from a transmitter located on a 
shore control station. A microwave system is commonly used 



19 



in this arrangement tut systems employing a laser can alsc 
be used for short-range work. Another LOP is generated by 
fixing an azimuth frcm a shore control station to the 
vessel. A second control station is used for an initial 
azimuth ty the observer. Azimuth determinations can be made 
after observing directions with a theodolite as an observer 
tracks the moving vessel. 

Ihere are twc ways to determine a range-azimuth 
position. The most common way is to have the theodolite and 
the transmitter occupy the same shore control station. 

Hence, the angle of irtersection, g, of the LOP's is always 
9Q0. This arrangement is commonly used by the National 
Ocean Service (NOS) fcr large-scale nearshore surveys. 

The other way is to have the theodolite and the 
transmitter occupy two different control points. Then the 
geometry is similar to that of the range-range position. 

The angle of intersection, g, is computed by trigonometric 
relationships among the azimuth of a line between the shore 
stations, the observed direction to the vessel, and the 
measured range to the vessel. 

4 . Azi muth - Azimuth 

Azimuth-azimuth positioning geometry is used for 
nearshore high-accuracy surveying. Theodolites are set over 
two control stations on shore. The vessel is sighted on 
simultaneously by the two theodolite observers, generating 
two visual LOP'S whose intersection define the vessel's 
location. Initial azimuths are fixed by sighting on control 
stations which are visible to the observers. 

The angle of intersection for an azimuth-azimuth 
position is dependent on the geometric relationships between 
the occupied stations, the initial stations, and vessel's 
position (Fig. 2.3). Assuming that theodolite observers 
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occupy stations 1 and 2 , and initial on stations 3 and 4, 



respectively, the observer at station 1 measures angle and 

the observer at station 2 measures y to the vessel. The 

2 

angle of intersection, 8 / is then computed by first deter- 
mining the forward arimuths, measured clockwise from the 
south, ficm stations 1 to 2 (a ), 1 to 3 (a ) , 2 to 1 

12 13 

( Cl ) , and 2 to 4 (a ) • The interior angles, 9 and 9 , of 

2 * 2 <* 12 

triangle 12P are 



9 - 

1 




+ Y. 




( 2 . 8 ) 



and 



9 => I a + Y - a I 

2 24 2 21 



(2.9) 



so the angle of intersection, 8 , at the vessel's location 
is 



8 = 180° - (9 + 9 ) 

1 2 



( 2 . 10 ) 



E. ClASEES OF EEBOBS 

All hydrographic positioning measurements are subject to 
error. The following sections discuss categories of errors 
and methods used to treat these errors. 

1 . Blunder s 

Blunders are gross mistakes which are generally due 
to the carelessness of the observer. Blunders can vary in 
magnitude, ranging from large errors which are easily 
detected, to small errors which may be barely distinguished. 
They can be detected by making repeated observations cr by 
carefully checking the data in the processing phase. 

Blunders occur in various forms and most can be avoided by 
carefully planning the data acquisition process. 
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Figure 2.3 Geonetry of an azimuth- Azimuth Position 

Consider the following as an example of a blunder 
associated with range-range geometry. An offshore pcsiticn 
is to be determined bj the intersection of two electronic 
LOP'S generated from transmitters located on known shore 
stations. The vessel is working west of a shoreline that 
runs generally in a north-south direction. As the hydrc- 
grapher faces the stations from sea, the southern shore 
station is mistakenly identified as left and the northern 
shore station as right. The resultant offshore position 
will plot to the east of the base line. This blunder is 
readily detected and can be easily remedied. 

Hot all types of blunders are so easily detected. 
Suppose an offshore position is to be determined by a 
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range-azimuth fix. fl range and an azimuth are generated 
from a known control station to the vessel's position. A 
second control staticn is used to fix the initial azimuth; a 
third shore control station is located 10 meters from the 
initial station and its coordinates are mistakenly used for 
the initial station in plotting. The resultant hydrographic 
position is in error, but this error will not be easily 
distinguished. 

Although most blunders have their origin ir human 
carelessness, some can be attributed to equipment malfunc- 
tion. For example, microwave systems which generate LOE's 
are known to become unsteady under certain conditions. 
Spurious range readings resulting from signal reflections 
can be recorded as true positioning data. In this case, the 
blunder may or may not be easily detected. 

In automated data acquisition systems, software has 
teen developed to detect the occurrence of anomalous range 
readings. By inputting a course and speed of a vessel trav- 
eling along a line, the computer can determine if the 
recorded position is valid based on the principle of dead 
reckoning. If the recorded position is found to be invalid 
the hydrographer will be immediately alerted to the situ- 
ation and can take action to remedy the problem. In non- 
automated systems the principle of dead reckoning is applied 
manually. Given the course and speed of the vessel, the 
validity of the position can be checked with spacing 
dividers. This involves checking the spacing between fixes 
recorded before and after the position in question. 

Before any type of error analysis is to he performed 
on the hydrographic positioning data, it is essential that 
all blunders be identified and properly treated. In 
general, careful planning coupled with thorough checking 
will minimize the occurrence of blunders. 
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Sy st ema t ic Errors 

Systematic errors occur with the same sign, usually 
of similar magnitude, and can te expressed in terms of a 
mathematical model. Systematic errors follow a defined 
pattern and occur in a number of consecutive related obser- 
vations. Repetition of measurements does nothing to mini- 
mize their effect. Ir the case of hydrographic positioning, 
systematic errors are identified and modeled by calibration 
of the measuring instrument against a known standard. Ihe 
following is a brief discussion concerning systematic errors 
and their treatment in relation to hydrographic positioning 
equipment. 

a. Theodolites 

In nearshore surveys the theodolite is used 
primarily for range-azimuth and azimuth-azimuth positioning. 
Systematic errors associated with the theodolite can be 
classified into two groups; those associated with the phys- 
ical design of the instrument and those involving the geom- 
etry of the positioning scheme. Some sources of systematic 
errors fEef. 8] associated with the physical characteristics 
of a theodolite are: 

i. The horizontal circle may be eccentric, 

ii. Graduations on the horizontal circle may not be 
uniform. 

iii. The horizontal axis of the telescope (about which 
it rotates) may not be perpendicular to the 
vertical axis of the instrument, 

iv. The longitudinal axis of the telescope may not he 
normal to the horizontal axis. 

V. The telescope axis and the axis of the leveling 
bubble may not be parallel. 
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These errors are usually small in magnitude and can be elim- 
inated by proper adjustment of the instrument by either the 
manufacturer or a qualified technician. 

The field hydrographer has ultimate control over 
the geometric systematic errors associated with a theodo- 
lite. In range-azimuth positioning the theodolite and 
transmitter may occupy the same horizontal control station. 
If the theodolite is not set directly over the station a 
resultant systematic error will occur in all measurements. 

It can be shown that these errors are non-linear but do 
follow a mathematical relationship. Likewise, if the trans- 
mitter is not located directly over the station, a similar 
type of lias occurs. Depending on the eccentricity of the 
theodolite, the vessel's range from the theodolite, and the 
scale of the survey — these errors can seriously affect the 
absolute accuracy of the offshore positions. 

In a similar fashion, it is also imperative to 
position the target directly over the horizontal control 
station used as an initial. Failure to do this will result 
in an error which will be propagated to offshore positions. 

Many situations arise in the field where it is 
advantageous to set a transmitter and theodolite over a 
single horizontal control station. Frequently it is 
feasible to construct a platform to accommodate both instru- 
ments; in a case where it is not, the position of an eccen- 
tric horizontal control station near the original station 
should be determined and that station used for the location 
of one of the instruments. The theodolite and the trans- 
mitter then occupy tvc known stations and the geometric 
source of systematic error is eliminated. 

b. Electronic Ranging Systems 

The systematic errors associated with electronic 
positioning systems are complex in nature and functions of 
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many variables. Munscn [Ref. 9, p. 4] addresses several 
problems associated with short-range systems used in hydro- 
graphic surveys. The most common problems with short-range 
systems are variation in range and calibration drift with 
time. Variations in internal eguipment time delays in the 
transmitter, the transponder, or the receiver can induce 
errors in measured ranges. For pulse systems such varia- 
tions can occur due to temperature dependence of components 
and fluctuations in signal strength at the transponder. 
Multipath effects are also a problem. Under some circum- 
stances a reflected wave and the directly transmitted wave 
arrive with a phase difference of ISO®. Cancellation or 
fading of the directly transmitted signal can result. 

NOS conducts base line calibrations of short- 
range positioning systems periodically during the course of 
a survey to minimize cr eliminate systematic error. In this 
process, a transmitter and receiver are each placed over 
control stations on shore and the measured range is compared 
to the true range. Ir this way the systematic error is 
eliminated by zeroing the instrument or by applying a 
constant correction to raw data. System checks are 
performed daily to assure there is no drift from the orig- 
inal calibration. A check can be accomplished by comparing 
a position defined by the ranging system to a known fixed- 
point position, to a sextant fix position, or an intersec- 
tion position. 

Munson [ Bef- 9, p. 5] also discusses sources of 
systematic errors associated with medium-range systems. Ihe 
most significant systematic errors occur as a function cf 
position due to varying propagation velocity. The medium- 
range electronic signal propagation velocity depends on the 
surface conductivity and transmission path (over water, over 
land, or over different types of land) . Because of this 
dependence, systematic errors as a function of position 
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occur at different effective phase velocities. Knowing the 
propagation velocity to use, or the phase correction to aake 
as a function of range, is a problem. Sky wave and storm 
interference also pose problems. At extreme ranges of cper- 
ation, sky wave interference can affect the more predictable 
ground wave, especially during nighttime operations. lane 
ambiguities are also a problem. Most systems are inherently 
ambiguous and must be zero set and continually mcnitcred for 
lane jumps or loss of signal which results in the loss cf 
lane count. 

NOS uses several techniques to determine the 
systematic error associated with medium-range positioning 
systems. These techniques involve determining a whole and 
partial lane count for phase comparison systems. Two of the 
more widely used techniques are comparison of three-point 
sextant fix positions to positions determined by the elec- 
tronic ranging system and calibration of the electronic 
system at a fixed point. In both techniques the whole lane 
counts are fixed by the calibration; correctors to the 
partial lane count are determined and applied to the raw 
ranging data. 

3 . Fa n dom Error s 

Eandom errors are chance errors, unpredictable in 
magnitude or sign, and are governed by the laws of prob- 
ability [Ref. 10, p. 1206], They are errors which remain 
after blunders and systematic errors have been removed. 

Random errors result from accidental and unknown combina- 
tions of causes and are beyond the control of the observer. 
Greenwalt [Ref. 12, p. 2] states they are characterized by: 

i. Variation in sign; positive errors occur with 
equal frequency as negative ones, 

ii. Small errors occur more frequently than large errors, 

iii. Extremely large errors rarely occur. 
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Eandom errors are unigue to specific types of pcsi- 
tioning eguipment and vary in magnitude depending on the 
precision of the instruments that are used. The following 
secticn outlines statistical methods for their treatment. 

C. IBIAiaEIlT OF RAHICfl EEBOBS. 

1 . Cne-Di me nsion al E rrors 

Certain basic statistical quantities must first be 
defined in the analysis of random errors. Consider a vessel 
moored securely to a fixed offshore platform. A number of 
ranges# n# from a microwave transmitter located on a shore 
control station are recorded. The mean of these observa- 
tions is 

n 

u = E — 

X i,l " (2.11) 

where x^ represents an individual observation. The standard 
error, s, of the observations is then 

^ - 1 ■ (2.12) 

where the quantity (x^- \i ) is referred to as the residual, 

1 X 

cr true error, v^ , of a particular observation. As n gets 
very large, the factor 1/n can be substituted for l/(n-1) in 
Equation 2.12. likewise, in treating the large sample, a 
can be substituted for s and m for , where y and g are 
the mean and standard error of the entire population. 

It is of interest to determine the probability of 
occurrence of a particular observation. The normal cr 
Gaussian distribution equation relates the residual of a 
particular random variable with the probability of its 
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1 



( 2 . 13 ) 



occurrence, and is given by 

P(v) = 






a \/^ 

Ibe plot of this equation yields the normal distribution 

curve (Fig. 2.4). The height of the curve above the 

vertical axis is proportional to the probability of a 

particular error occurring'. 

The probability of a residual falling between any 

two residuals v and v can be computed by integrating 
1 2 

Equation 2.13 as 

.2 



P(V) = 



1 

\fZTi 






dv 



(2. 14) 



P(V) 




I 

! 






Figure 2.4 The Hormal Distribution 

This integral is difficult to evaluate analytically 

so tables have been compiled to aid in computations. Fcr 

V = +0 and v = -a , it can be shown that P(v) = 0.6827. In 
1 2 

ether words, the protability that a particular observation 
will fall within + la of the mean is 68.27 percent- 
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Beturning to the example of the vessel moored tc the 
offshore platform, the mean and the standard error for the 
observations are easily computed. With this information and 
Equation 2.14, the probability of a range error falling 
within specified limits can be computed. Conversely, by 
fixing a probability, the associated limits of the range 
error can be computed . In statistical terms, a particular 
observation will fall within specified limits with a certain 
conf i dence . 

Actual values of one-dimensional standard errors for 
hydrographic positioning equipment are a subject of debate 
between manufacturers and users. Some manufacturers of 
microwave positioning equipment claim standard errors of +1 
meter. Cn the other hand, Munson [Ref. 9, p. 6] states that 
microwave systems demonstrate accuracies of 3 meters at 
short ranges but show larger errors at ranges of 15 km and 
greater. NOS assumes a 3-meter standard error in all of its 
short-range accuracy ccmputaticns. It is apparent that 
further study is needed to adequately define the nature of 
errors associated with electronic positioning equipment. 

Waltz [Ref. 13] performed an extensive study to 
determine the pointing error of a Wild T-2 theodolite. His 
results showed that the pointing error associated with this 
Instrument under hydrographic survey conditions was about 
1.3 meters and was independent of distance. 

2 . I wo-Di mensi onal Err ors 

The intent of this paper is to apply statistical 
methods developed by ethers tc a hydrographic data set 
containing two-dimensional errors which are defined by two 
randem variables. Lengthly and complex derivations are not 
presented. Burt [Hef. 3] and Heinzen [Ref. 2] show adequate 
derivations of formulas associated with two-dimensional 
errors and can be referenced for full details. 
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Ihe followicg assumptions are made concerning tvio- 
dimensicnal errors associated with intersecting LOP's: 



i. Ihe random errors of each LOP are normally 
distributed. 

ii. Systematic or bias errors have been removed from 
the observations. 

iii. The intersecting LOP's are coplanar. 

iv. The error LOP’s are parallel to the exact LOP's. 



In developing a usable mathematical model for accuracy 
determinations, the four assumptions hold to a high degree 
for all hydrographic positioning geometries. 

Consider again the vessel moored to a fixed offshore 
platform. Assume two ranges are measured from two different 
shore control stations at the same time and that the range 
readings are uncorrelated. The observation of this pair of 
ranges is repeated many times. After a large number of 
observations, the means and standard errors of the indi- 
vidual ranges are determined. Suppose the mean ranges, or 
the actual LOP's, intersect at an angle of 90 ° and that the 

computed standard errors are egual (a =a )• If each data’ 

1 2 

pair plotted, the spread of points about the mean 

coordinates results in a circular cluster (Fig. 2.5). A 
higher density of points occurs near the intersection of the 
mean ranges and the density of points decreases outward from 
the intersection of the mean ranges. 

In this special case, which is called a circular 
normal distribution, the probability of a point falling 
within a specified radius, R, from the intersection of the 
mean ranges is 



P(R) » 1 - e 






3 1 



(2. 15) 



where a = o = is deiiced as the circular standard 

1 2 ^ 

error. Dsing Eguaticr 2. 15, R can be computed by fixing 

P (R) , cr ccnversly, P (R) can be computed by fixing R. 

Letting R = o = a = a / then P(R) = 0.3935. In other 
1 2 c 

words, 39.35 percent cf all errors in a circular normal 
distribution are not expected to exceed the circular stan- 
dard error [Ref. 12, pp. 25 t26]. 




Figure 2.5 Circular Normal Distribution 

In the case where the two uncorrelated LCP’s inter- 
sect at an angle other than 9Qo or a , the contours cf 

I 2 
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equal density are ellipses centered about the point defined 
by the intersecting ICP's (Fig. 2.6). The two-dimensional 
probability density function becomes [Ref. 1, p. 136] 






K" 

7 



(2. 16) 



LOP 2 • 




Figure 2.6 Error Ellipse Formed by Two Oncorrelated ICP's 
where 

V is the residual in the direction of the semi-major axis 

X 

of the error ellipse/ 

Vy is the residual in the direction of the semi-minoi 
axis. 
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is the standard error in the direction of the semi- 
major axis, 

a is the standard error in the direction of the semi- 

y 

miner axis. 



and 




(2. 17) 



The solution of Equation 2.16 with values of K for different 

P's yields the results in Table I [Bef. 12, p. 23]. For a 

39.35 percent probability, the axes of the ellipse are 

1.0000a and 1.0000 a ; for a 50 percent probability, the 

axes are 1. 1774 a an<J^ 1. 1774 cr • 

X y 



TABLE 


I 






Values of the 


Constant K 






PROEAEILITY 


K 






39.35% 


1.0000 






50.00% 1 


1 . 1774 






63.21% 


1.4142 






90.00% 


2.1460 






99.00% 


3.0349 






99.78% 


3.5000 




1 



The error ellipse can be used for accuracy computa- 
tions by developing relationships for a and a in terms of 

X y 

the initial information a / a , and g . Bowditch [Bef. 10, 

l 2 

p. 1213] gives the following equations for independent lOP's 
relating these quantities: 



= 



2?!^ ( *■ t -J (cr; + cM‘ - 4 sln ‘6 ) ,2. 18) 



34 



and 




Id these equations, Q is assumed to be the acute angle 
betwesD the LOP’s. 



on more manageable fcrms. In range-range and azimuth- 

azimuth positioning it is often assumed that a = a = a. 

1 2 

Iguations 2.18 and 2-19 then reduce to 



In the concentric range-azimuth case, c * a , and S 

1 2 

eguals 90°. Equations 2.18 and 2.19 then simplify to 



where o > o and a ^ a • 

12 X y 

The case for correlated LOP’s is more complex. The 



tion from a linear shewed coordinate system to an uncorre- 
lated rectargular coordinate system. The following 
discussion is taJcen from Heinzen [Ref. 2, pp. 49-53]. 

Assume a hydrographic position is established by the 
intersection of two correlated LOP’s (Fig. 2.7a). LOP 1 and 



In certain special cases, the above equations take 



“ 2sin(>56) ° 



( 2 . 20 ) 



and 




( 2 . 21 ) 




( 2 . 22 ) 



and 




(2.23) 
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A 
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RECTANGULAR 

B 



UNCORRELATED 

rectangular 

C 



Figure 2.7 Coordinate Transformations for Correlated IOP*s 



10? 2 are the coordinate axes in the skewed coordinate 

systeB, with standard errors a and a . The semi-major and 

1 2 

semi-minor axes of tie error ellipse are not coincident with 
the skewed coordinate system axes. The correlation coeffi- 
cient between the two LOP's is p . Assume (j > a • 

12 12 

The standard errors and correlation coefficient in a 
correlated rectangular coordinate system with axes A and B 
must now be determined. A coordinate transformation from 
the skewed system to the correlated rectangular system must 
be made yielding the standard errors along the new coordi- 
nate axes (Fig. 2.7b) 



sin'^B 



(a^ + 2p a a cosB + o^) - 
1 12 1 2 2 



(2.24) 
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and 



Ob = 02 



(2.25) 



Ihe correlation coefficient in the correlated rectangular 
system is 

0 0 0 2 
) { 1 + D 7-t^)cOS 6 + (-r^) 

1 2 



°ab ' ' 1 ^ (5^)cos=6>-'^ 

1 1 1 ' ' 



lo determine o and 0 , a second coordinate transformation 
X y 

must be performed from the correlated rectangular system to 
an uncorrelated rectangular system with axes X and I (Fig. 
2.7c). Ihe semi-majoi and semi-minor axes of the error 
ellipse are then 



^ * 






1 + a 1 - 









T^TT^ 



(2.27) 



and 






(2. 28) 



When 0,. = 0, these equations become identical to the 
simplified versions in Bowditch £Ref. 10]. 

Ihe orientation of the semi-major and semi-miner 
axes relative to the intersecting LOP's is the third param- 
eter which fixes the error ellipse. The angle 0 (Figs. 2.6 
and 2.7) is measured counter-clockwise from LOP 1 to the 
semi-major axis of the error ellipse [Ref. 11] and is given 
by 



, o^sin( 2 B) + 2 p 00 s 1 n( 8 ) 

9 = i arctan { — r 

^ 0^cos(28) + 2 p 00 cos(S) + 0 

1 12 i 2 2 



} 



(2.29) 



For the special case cf 



0 =0 and 

1 2 



P =0, 
12 



8 

I 



(2.30) 
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The orientation of the error ellipse in an orthogonal coor- 
dinate system can be represented by adding or subtracting 9 
to the orientation of LOP 1 . Care must be taken on deter- 
mining the guadrant of the outcome. As a general rule, the 
error ellipse always lies within the acute angles formed by 
the intersecting LOP’S. 

The orientation and dimensions of the error ellipse 
provide a useful index for evaluating the accuracy of a 
hydrographic position. Its greatest attribute is that it 
accurately represents the error distribution about the 
intersection of two TCP’s in terms of a fixed probability. 

It is interesting to examine the variation in the relative 
dimensions and orientations of error ellipses as they vary 
in a range-range configuration with o = a = o (fig- 2.8) . 

I 2 

The dimensions of the ellipses are specified by Equations 
2.20 and 2.21 and o and a are functions of g only for 

X y 

fixed o . Therefore, the dimensions of the ellipses remain 
constant along a contour of constant g ; only the 
orientation changes. A line of constant 6 is a circle 
which includes stations L and B. Note that the dimensions 
of the ellipses for S ’ s of 30 ° and 150 ° are identical. The 
ellipses about the 90 ° angle of intersection contour are 
circles and represent the strongest possible positions in 
this scheme. With varying 8’s, the directional nature of 
the distribution can he noted. 

~ • Cir cul ar P re c ision I ndexes 

Although the error ellipse gives a true representa- 
tion of the error distribution about a hydrographic posi- 
tion, its use has certain drawbacks. The characteristics of 
the ellipse must be specified by the three quantities a , 

X 

o^, and 9. A single figure for evaluating the positional 
accuracy cannot be used. Greenwalt [fief. 12 , p. 26 ] states 
that when 0 and q are not equal, a circular error 

X V 
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1 




Figure 2.8 Error Ellipses around a Range-Range System 
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distributicn can be substituted for the elliptical distribu- 
tion. This substitution can be satisfactory for error 
analysis vithin certain ratios. However, when this 

ratio is small the distortion introduced by the circular 
distribution may beccie misleading. 

a. Root Mean Square Error 



The terms radial error, root mean square error, 

and d are identical in meaning when applied to two- 

rms 

dimensional errors [Ref. 10, p. 1229 ]- The term d is 
defined as the square root of the sum of the squares of the 
standard errors along the major and minor axes of the error 
ellipse. That is 



d 



= \la^ + 
rms '' X 



(2.31) 



where and o are given by Equations 2.18 and 2.19. 

A more direct form of 2.31 is given by [Ref- 2, p. 54] 



d 



rms 



1 

sinS 





(2.32) 



for uncorrelated LOP’s. For range-range and azimuth-azimuth 

positioning, with a = o = a , Equation 2.32 reduces to 

1 2 

‘^rms " (2.33) 

For range-azimuth positioning, 8 = 90° and Equation 2.32 
becomes 



d 



rms 



\Ja^ + 0^ 

l 2 



(2.34) 



The more general form of Equation 2.32 for both correlated 
and uncorrelated LOP’s [Ref- 2, p- 59] is 



nns 



sine 






+ + 2 p 00 cosS 

2 12 1 2 



(2.35) 



where P is the correlation coefficient. 
12 
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An error circle with a radius of one d can be 

rms 

constructed about the intersecting LOP’s (Fig. 2.9). Two 

d is the radius cf the error circle obtained using two 
ms 

times the values of a and a in Equation 2,31. For an 

X y 

elliptical error distribution, the probability associated 

with a specific value of d varies as a function of the 

ms 

eccentricity of the error ellipse (Table II) . The prob- 
ability associated with one d varies from 63-2 percent to 
68.3 percent, while the probability associated with two d 

nns 

varies between 95.4 percent and 98.2 percent. 







Figure 2,9 The d Error Circle 

rms 



NOS uses d^„ 
rms 

Umbach [Eef, 14, p- 4-25] 
direct distance measuring 
the value cf d^^ is less 



as an accuracy specification, 
states that super high frequency 
systems would be used only when 
than or equal to; 
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i • 


0.5 mm at the scale of 
1:20,000 and smaller. 


the 


survey 


for 


scales of 


ii . 


1.0 mm at the scale of 
scale surveys, or 


the 


survey 


for 


1:10, 000 


iii. 


1.5 mm at the scale of 
1:5,000 and larger. 


the 


survey 


for 


scales of 




The major advantage 


of 


using d 


rms 


as a precision 



index is its ease of computation. Some hydrographers draw 

analogy tetween the varying probability associated with one 

d (63.2 percent to 68.3 percent) and the fixed prcb- 
rms 

ability associated with a one-dimensional standard error 

(68.3 percent). In fact, d has very little statistical 

nas 

meaning. Ihe obvious problem with using as a precision 

index is the varying probability associated with the error 
circle. For this reason Greenwalt [Ref. 12, p. 31] recom- 
mends against its use. 



TftBlE II 

Probabilities Associated With d 

ms 

PROBABILITY 



a 

y 


a 

X 


LENGTH OF 


f\ 

cn 


0.0 


1.0 


1. 000 1 


0.683 


0. 1 


1.0 


1.005 


0.682 


0.2 


1.0 


1.020 1 


0.682 1 


0.3 j 


1.0 


1.042 


0.676 


0.4 


1.0 


1.077 1 


0.671 


0.5 


1.0 


1. 118 


0. 662 


0.6 


1.0 1 


1. 166 1 


0.650 


0.7 


1.0 1 


1. 220 


0.641 


0.8 1 


1.0 


1.280 i 


0.635 


0.9 


1.0 


1. 345 


0. 632 


1.0 ] 


1.0 


1.414 J 


0.632 



2 d 



rms 



0. 954 
0. 955 
0.9 57 
0. 961 
0. 966 
0. 969 
0. 973 
0. 977 
0. 980 
0. 981 
0.982 
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b. Circles cf Eguivalent Probability 



Burt £Sef. 3] presents a method for translating 
ellipses of equivalent probability into circles of equiva- 
lent probability. Tc utilize this method, it is first 
necessary to compute the eccentricity of the error ellipse, 
c, by the equation 



a 




(2.36) 



where a > a . 

X y 

Harter [ Eef . 15] compiled Tables III and IV 
which are taken from Bowditch [Eef. 10, p. 1215]- Harter's 
data are given in terns of the eccentricity, c, a parameter, 
K, and a probability, P. The parameter, K, when multiplied 
by gives the value of the radius, R, of the circle of the 
corresponding probability shown in Table III. That is, 

(2.37) 

The probability of a point falling inside a 
circle cf specified radius can be computed by entering Table 
III with c and K as arguments. Given a fixed probability, K 
is determined by entering Table IV using c and P as argu- 
ments. The radius of the probability circle is then 
computed using Equation 2.37. 

Using confidence ellipses has certain advantages 
over confidence circles of equal probability. First, the 
directicnal nature of the true error distribution is net 
represented in the confidence circle method even though both 
methods give an accurate measure of confidence. Second, the 
area of the confidence ellipse is always less than or equal 
to the area of the confidence circle. The area of a 
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TABLE III 

Probabilities, Giyen c and H 



c 

K \ 


a 0 


0. 1 


0. 2 


a 3 


0. 4 


a 5 


0. 6 


0. 7 


0. 8 


0. 9 


1. 0 


0. 1 


. 0790557 


. 0443987 


. 0242119 


. 016417G 


. 0123875 


. 0099377 


. 0082940 


. 0071157 


. 0062299 


. 0055400 


. 0049875 


0. 2 


. 1585194 


. 1339783 


. 0884533 


. 0628396 


.0482413 


.0390193 


. 0327123 


. 0231415 


. 0246824 


. 0219757 


. 0198013 


0. 3 


. 235S228 


. 2213804 


. 1739300 


. 1318281 


. 1039193 


. 0851535 


. 0719102 


. 0621386 


. 054G598 


. 04S7639 


. 0440025 


0. 4 


. 3108435 


. 3010228 


. 2635181 


. 2139084 


. 1742045 


. 1451808 


. 1237982 


1 . 107G237 


. 0950495 


. 0S.5032G 


. 0768837 


0. 5 


. 3829249 


. 3755SS4 


. 3481790 


. 3003001 


. 2532953 


. 2152880 


. 1357448 


. 1020829 


. 1443941 


. 1290286 


. 1175031 


0. 6 


. 4514938 


. 4457708 


. 4255605 


. 3S46374 


. 33573S4 


. 2914682 


. 2548177 


.2251114 


. 2009797 


. 1811783 


1647298 


0. 7 


. 5160727 


. 5115048 


. 49G0683 


. 4633258 


. 4170862 


. 3699305 


. 3280302 


. 2925654 


. 2629373 


. 2.381.583 


. 2172955 


0. 8 


. 5762892 


. 572.5957 


. 5604457 


, 5349387 


. 4941882 


. 4474207 


. 402.5628 


. 3627122 


. 32S3453 


. 2989700 


1 . 2738.510 


0. 9 


. 6318797 


. 6288721 


. G191354 


. 5993140 


. 5651564 


. 5213998 


. 4759375 


. 4333628 


. 3953279 


. 36.20135 


. 3330232 


1. 0 


. 6826895 


. 6802325 


. G723580 


. 6568242 


. 6291249 


. 5900953 


. 5461319 


. 5025790 


. 4621421 


. 4257553 


. 3934693 


1. 1 


. 7286679 


. 7266597 


. 7202682 


. 7079681 


. 6859367 


. 6524489 


. 6116316 


. 5687467 


. 5272462 


. 48S7S73 


. 4539256 


1. 2 


. 7698607 


. 76S2215 


. 7630305 


. 7532175 


. 7359558 


. 7079973 


. 6714269 


. 6306168 


. 5893494 


. 5498736 


. 5132177 


1. 3 


. 8063990 


. 8050G48 


. 8008554 


. 7929968 


. 7793550 


. 7567265 


. 7249673 


. 6873122 


. 6474394 


. 6079822 


. 5704426 


1. 4 


. 83S4867 


. 8374049 


. 834001S 


. 8277048 


. 8169851 


. 7989288 


. 7720S89 


. 7383089 


. 7007900 


. 6623035 


. 6216889 


1. 5 


. 8GG3856 


. 8655127 


. 8G27728 


. 8577362 


. 8493071 


. 8350816 


. 8129287 


. 7833962 


. 7489500 


. 7122546 


. 6753473 


1. 6 


. 8904014 


. 8897008 


. 8875060 


. 8834914 


. 8768644 


. 86575.59 


. 8478393 


. 8226246 


. 7917194 


. 7574708 


. 7219627 


1. 7 


. 9108691 


. 9103102 


. 9085619 


. 90537G6 


. 9001746 


. 8915536 


. 8773116 


j . 8.562471 


. 8291137 


. 7977882 


7462539 


1. 8 


. 9281394 


92769G4 


. 92G3125 


. 9237989 


. 9197275 i 


< . 9130680 


. 9019110 


•. 3346624 


. SO 13238 


. 8332175 


. 5021013 


1. 9 


. 94256G9 


. 9422182 


. 9411299 


1 . 9391586 


. 9359855 


i . 9308615 


. 9222277 


. 90S3609 


. 8.886731 


. 8639149 


. 8355255 


2. 0 


. 9544997 


. 9542272 


. 9533775 


- 9518415 


. 9493815 


. 9454546 


. 9388418 


. 9278799 


. 9115762 


. 8901495 


. 8G46647 


2. 1 


. 9G42712 


. 9640598 


. 9634011 


. 9622127 


. 9603170 


. 9573205 


. 9522999 


. 9437668 


. 9305013 


. 9122714 


. 8.H97495 


2. 2 


. 9721931 


. 9720304 


1 . 9715237 


. 9706109 


. 9091597 


. 9608845 


. 9631017 


1 . 956.5522 


. 9459386 


I . 9306821 


. 9110784 


2. 3 


. 9785518 


. 9734275 


. 9780408 


. 9773450 


1 . 9762419 


. 9745239 


. 9716934 


9667306 


. 9583739 


. 9458085 


. 928994G 


2. 4 


. 9S36049 


. 9835108 


9832180 


. 9826918 


. 9S18594 


. 9805703 


. 9784651 


.9747495 


.9682698 


1 . 9.580804 


. 942S652 


2. 5 


. 9875807 


. 9875100 


. 9872900 


. 9868953 


. 9862720 


. 9853112 


. 9837569 i 


. 9810035 


.9760522 

1 


1 . 9679136 


I . 9SG0G31 


2. 6 


. 9906776 ' 


. 9906249 


. 9904612 


. 9901674 


. 9897045 


. 9889934 


. 9878.527 


. 98.58331 


. 9821023 1 


1 . 97.56959 


1 . 96.5952.5 


2.7 


9930661 


. 9930271 


. 9929062 


. 9926894 


. 9923483 


. 9918260 


. 9909944 


. 9S95268 


. 9867530 I 


1 . 9817837 


! . 97387S6 


2. 8 


. 9948S97 


. 9948612 


. 9947727 


. 9946141 


. 9943649 


. 9939842 


. 9933821 


. 9923249 


. 9902888 


. S864876 


. 9o01589 


2. 9 


. 9962684 


. 9962477 


. 9961834 


. 9960684 


. 9958878 


. 9956126 


. 99.51798 


. 9944246 


.99294.S2 I 


! . 99CXJ803 


. 9 <5.5 0792 


3.0 


. 9973002 


. 9972853 


. 9972391 


. 9971564 


. 9970266 


. 9968294 


. 9965205 


. 9959854 


.9949274 1 


[ . 9927925 


. 9888910 


3. 1 


- 9980648 


. 9980542 


. 9980212 


. 9979622 


. 9978699 


. 9977296 


. 9975109 


. 9971348 


. 99638.51 


. 9948168 


. 991S113 


3. 2 


. 9986257 


. 9986182 


. 9985949 


. 9985533 


. 9984880 


. 9983892 


. 99823.56 


. 9979733 


. 9974478 I 


; . 9963105 


. 9940240 


3. 3 


. 9990332 


. 9990279 


. 9990116 


. 9989824 


. 9989368 


. 9988677 


. 9987607 


. 9985792 


.9982147 ! 


' . 9974004 


. 9956822 


a 4 


. 9993261 


. 9993225 


. 9993112 


. 9992909 


. 9992.593 


. 9992115 


. 999 1376 


. 9990129 


.9987626 


1 . 9981S68 


. 9969113 


3. 5 


. 9995347 


. 9995323 


. 9095245 


. 9995105 


. 9994888 


. 9994559 


. 9994053 


. 9993204 


. 9991502 I 


. 99S7480 


. 9978125 


3. 6 


. 9996818“ 


. 9996801 


. 9996748 


. 9996653 


. 9996505 


. 9996281 


. 9995938 


. 9995364 


. 9994218 


. 0991442 


. 9984662 


3. 7 


. 9997844 


. 9997832 


. 9997797 


. 9997733 


. 9997633 


. 9997482 


. 9997251 


. 9996867 


.9990102 


. 9994208 


. 9989352 


3. 8 


. 9998553 


. 9998545 


. 999 So 22 


- . 9998478 


. 9908412 


. 9998311 


. 99981.57 ' 


. 9997902 


. 9997396 1 


1 . 9996119 


. 9992G82 


a 9 , 


. 9999038 


. 9999033 


. 9999018 


. 9998989 


. 9998945 


. 9998S78 


. 9998776 1 


. 999860S 


.9998276 ! 


. 9997426 


. 9995020 


4 . 0 1 


. 9999367 


. 9999363 


. 9999353 


. 9999334 


. 9999305 


. 9999261 


. 9999195 


. 9999085 


. 0998S70 


. 9999309 


. 999664.> 


4. 1 


. 9999587 


. 9999585 


. 9999578 


. 9999566 


. 9999547 


. 9999519 


. 9099475 


. 9999404 


. 9999266 


. 9998900 * 


. 9997763 


4 . 2 


. 9999733 


. 9999732 


. 9999727 


. 9999720 


. 9999707 


. 9999689 


. 9999661 


. 9999616 


. 9999527 


. 9999292 


. 9998523 


4. 3 


. 9999829 


. 9999828 


. 9999826 


. 9999821 


. 9999813 


. 9999801 


.9999753 ' 


. 9999754 


. 9999098 


. 9999.548 


. 9999034 


4 . 4 


. 9099892 


. 9999891 


. 9999889 


. 9999S86 


.9999881 


. 9999874 


. 9999863 


. 9999845 


. 9999809 


. 9999715 


. 999937.5 


4. 5 


. 9999932 


. 9999932 


. 9999931 


. 9999929 


. 9999925 


. 9999921 


. 9999914 


. 9999902 


. 9999881 


. 9999822 


. 9909599 


4 . 6 


. 9999058 


. 9999957 


. 9999957 


. 9999955 < 


. 9999954 


. 9999951 


. 9999947 


. 9999939 


. 9999926 


. 9999S89 


. 9999746 


4 . 7 


. 9999974 


. 9999974 


. 9999973 


. 9999973 


. 9999971 


. 9999970 


. 9999967 


. 9999963 


. 99999.5.5 


. 9999932 


. 9999'^40 


4. 8 


. 9999984 


. 9999984 


. 9999984 


. 9999983 1 


. 9999983 


. 9999982 


. 99999 SO 


. 9999977 i 


. 9999972 


. 9999959 


. 9999901 


4 . 9 


. 9999990 


. 9999990 


. 9999990 


. 9999990 


. 9999990 


. 9999989 


. 9999988 


.9999986 


. 99999S3 


. 9999975 


. 9999939 


5. 0 


. 9999994 


. 9999994 


. 9999994 


. 9999994 


. 9999994 


. 9990993 


. 9999993 


. 9999992 | 


. 9999990 


. 99999S5 


. 9999963 


5. 1 


. 9999997 


. 9999997 


. 9999997 


. 9999906 


. 9999996 


. 9999996 


. 9999996 


. 9999995 ! 


. 9999994 


. 9999991 


. 9999978 


5. 2 


. 9999998 


. 9999998 


. 9999998 


. 9999998 


. 9999998 


. 9999998 


. 9999998 


. 9999997 I 


. 9999997 


. 9999995 


. 9999987 


5. 3 


. 9999999 


. 9999999 


. 9999999 


. 9999999 


. 9909999 


. 9999999 


. 9999999 


. 9999998 


. 9999998 


. 9999997 


. 9999992 


5. 4 


. 9999999 


. 9999999 


. 9999999 


. 9999099 


. 9999999 


. 9999999 


. 9999999 


.9999999 


. 9999999 


. 9999998 


. 9999995 


5. 5 


1. OOOOOOO 


1. 0000000 


1. 0000000 


1. 0000000 


1. 0000000 


1. 0000000 


1. 0000000 


. 9999999 


. 9999999 


. 9999999 


. 9999997 


5. 6 






i 










1. 0000000 


1. OOOOOOO 


. 9999999 


. 9999998 


5. 7 




















1. OOOOOOO 


. 9999999 


3. 8 

5. 9 

6. 0 






1 




1 












1. OOOOOOO 
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TABLE 17 

Radii of Circles Given c and P 



c 

p 


0. 0 


1 

a 1 

1 


1 

a2 


a 3 


0. 4 


0. 5 


0. 6 


0. 7 


0. 8 


0. 9 1 


1.0 


.5000 

.7500 

.9000 


a 67+49 
1. 15035 
1. 644S5 


0. 68199 

1. 15473 
1. 64791 


0. 70585 

1. 16825 
1, 65731 


0. 74993 

1. 19246 
1. 67383 


0. 80785 

1. 23100 
1. 69918 


0. 87042 

1. 28534 
1. 73708 


0. 93365 

1. 35143 
1. 79152 


0. 99621 

1. 42471 

1. 8G253 1 


1. 05769 1 
1. 50231 
1 1. 94761 

1 


i 1.11807 
1.58271 
2. 04236 


1. 17741 
1.06511 

2. 14597 


.9500 

.9750 

.9900 


1. 95996 

2. 24140 
Z 57583 


1. 96253 

2. 24365 
2. 57778 


1. 97041 

2. 25053 
2. 58377 


1. 98420 

2. 26255 
2. 59421 


2.00514 
2. 28073 ; 
2. 60995 


2. 03586 
2. 30707 
2.63257 


2.08130 
2. 34581 
2. 66533 


2. 14598 
2.40356 
2. 71515 


2.23029 
1 2.48494 
2. 79069 


2. 33180 1 
, 2,58999 i 
1 2.89743 


: 2. 44775 
, 2.71620 
1 3.03485 


.9950 

.9975 

.9990 


2. 80703 ' 

3. 02334 
3.29053 ' 


' 2. 808S3 
Z 02500 
' 3. 29206 


2. 81432 

3. 03010 
3. 29673 


2. 83289 

3. 03898 
3. 30489 


2.83830 ; 
3. 05234 ! 
3.31715 


j 2.858-94 ' 
1 3.07144 
3. 33464 


2. £8859 

3. 09871 
3. 35949 


2. 93347 

3. 13969 
3. 39647 


3. 00431 
3. 20556 
3. 45698 


1 3. 1 1073 
1 3. 31099 1 
! 3.55939 1 

1 1 


3. 25525 
1 3. 46164 
1 3. 71692 



1 



confidence ellipse is 






(2.38) 



where K is the appropriate prctability conversion factor 
(Table I) . The area of the SO percent confidence circle is 



A 



c 



= itR2 



(2.39) 



where B is given by Eguation 2.37. For a condition where 

0=0=3 meters, and S = 30°, the area of the 90 percent 

1 2 

confidence ellipse is 261 square meters, while the area ox 
the confidence circle is 587 square meters. For both stan- 
dard errors equaling 10 meters and 6 = 30°, the 90 percent 
confidence ellipse has an area of 921 square meters and the 
confidence circle has an area of 2894 square meters. From 
an operational perspective, the difference in areas between 
ellipses and circles have significant implications which 
will be discussed in Chapter v. 

The following examples are presented to deiicn- 
strate methods fcr ccmputing the parameters of error 
ellipses ard confidence circles for several hydrographic 
positioning geometries. 
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E xa m pie J 

A vessel is conducting a hydrographic survey 

using range-range geometry. The two LOP's generated ty 

microwave transmitters have standard errors of a = 3 meters 

1 

and o =4 meters. The angle of intersection 6 at the 
2 

vessel is 30°. Assume the LOP’s are un correlated . Compute 
the probability that the vessel's position will be within a 
circle of 10-meter radius with the center at the intersec- 
tion of the LOP’S. 

Recalling Eguations 2.18 and 2.19, the values of 

o and o are found to be 9.79 meters and 6.14 meters, 

X y 

respectively. From Eguation 2.36 

a 

c = = 0.633 

o 

X 

and frcffi Eguation 2.37, with R = 10 meters, 

K = 1.032 

Entering Table III and using interpolated values for c and 
K, the probability that the vessel’s position will be within 
a circle of 10-meter radius centered at the intersection of 
the ICP’s is 



P = 53.2% 



Examp le 2 

A vessel is conducting a hydrographic survey 
using range-azimuth cecmetry. The range LOP generated hy 
the microwave transmitter has a standard error of 3 meters. 
The azimuth LOP determined by theodolite observation has a 
standard error of 1.3 meters at all ranges. Compute the 
radius of the 90 percent confiderce circle at the vessel's 
position. 

In the range-azimuth case 0 = 90° and the ICP's 
are unccrrelated . Therefore, 

o =0 =3.0 meters 

1 X 

and 
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a 



1. 3 meters 



then 



a 

2 



c = -^ = 0.433 

Table IV is entered with the values of ? = 0.9 and c = 
0.433 . The value for K is found to be 



K = 1.7117 

Using Equation 2.37, the radius of the 90 percent prob- 
ability circle is fourd to be 

R = 5. 14 meters 

The probability that the vessel's position will be within a 
circle of 5.14-meter radius centered at the intersection of 
the ICP's is 90 percent. 



Example 3 



A vessel is conducting a hydrographic survey 
using hyperbolic-hyperbolic gecmetry. The hyperbolic LC? 
generated by the 1.6-t?Hz electronic positioning system has a 
standari error of 0.05- lane on the base line. The correla- 
tion coefficient { p ) between the two LOP's is known to be 

1 2 

0.4. Compute the radius of the 90 percent confidence circle 
at the vessel’s position. 

The rectangular plane coordinates of the master 
(M) , two slaves (G and R) , and the vessel's position (P) are 

I COORDINATE Y COORDINATE 



(m) 



(m) 



R 172,679.1 

G 308, 679.1 

M 241,738.2 

P 223,172.5 



62,540. 4 

98.540. 4 

21.325. 4 
169, 264. 2 
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Given the frequency cf 1.6 MHz, X = 187.37 meters frcm 

Equation 2.3. The lare width along the base line is 

= 93.68 meters from Equation 2.5. Using the law of cosines 

from plane geometry, the subtended angles a and a are 

32.470 and 43.25°, respectively. The angle of intersection 

of the two hyperbolas at P is 37.86° from Equation 2.7. The 

lane widths at P are w = 254. 19 meters and w = 335. 06 

^ g 

meters from Equation 2.6. The standard errors of the green 

(0 ) and red (a ) hyperbolas, respectively are a = w g = 
12 1 g base 

16.7 meters and a = w a = 12.7 meters. These standard 

2 r base 

errors are in a linear skewed coordinate system and must be 

transformed to an unccrrelated rectangular system. Frcn 

Equations 2.18 and 2.19, the values of a and o are 36.9 

a b 

meters and 12.7 meters, respectively. The correlation coef- 
ficient in the correlated rectangular system ) is then 

0.737 frcm Equation 2.26. The semi-major and semi-minot 
axes in the uncorrelated rectangular system are 38.1 meters 
and 8.3 meters, respectively, from Equations 2.27 and 2.26. 
The eccentricity is 

a 

c = ^ = 0.218 

X 

Table IV is entered with the values of P = 0.9 and c = 

0.218. The value for K is found to be 

K = 1.6602 

From Equation 2.37, the radius of the 90 percent probability 
circle is found to be 

R = 63.3 meters 

The probability that the vessel's position will be within a 
circle of 63.3-meter radius centered at the intersection cf 
the LCP's is 90 percent. 
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III. EXPERIBENT MSIGN AND IMPLEMENT ATIO N 



The goals of this chapter are to demonstrate that hydro- 
graphic positioning accuracy can be classified based or the 
radii cf 90 percent confidence circles determined by using 
Burt’s method and to show that, based on the same criteria, 
accuracy predictions can be made for survey planning 
purposes. 

A. BATA ACQUISITION IBOCEDUBES 

The data used for analysis and prediction consisted of 
range-range, azimuth-azimuth and range-azimuth survey infor- 
mation. The data were acquired by Naval Postgraduate School 
<NPS) students in a Hydrographic Sciences course. Although 
the course was structured as a training exercise, the data 
acguisiticn procedures utilized were nearly identical to 
those which are practiced by NOS. 

A total of 453 hydrographic positions were recorded 
during the survey of a nearshore area in southern Monterey 
Bay, California. Of the positions used for analysis, 2S2 
were range-range, 81 were range-azimuth, and 80 were 
azimuth-azimuth. All survey information was recorded by 
hand in sounding volumes. The vessel used was a 36-foot 
Uniflite with a fiberglass hull and twin engines. The 
survey was conducted on October 28, November 16, 23, and 30, 
1983. Electronic ccntrol and calibration stations used for 
the survey included C£E MON 1978, MUSSEL 1932, BEACH LAB 
1982, MCKTEF.EY AMERICAN CAN COMPANY STACK 1932, MONTEREY 
RADIO STATION KMBY MAST 1962, MONTEREY HARBOR LIGHT 6 1978, 
and MCNTEEEY BLUE LIGHTHOUSE (Fig. 3.1). With the exception 
of MONTEREY BLUE LIGHTHOUSE, which is a low-order position. 
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Figure 3. 1 Hydrographic Survey Area 

all stations are of third-order or better and are published 
in the National Geodetic Survey Data Base. 

For azimuth-azimuth and range-azimuth positioning, 
azimuths were measured with a Wild T-2 theodolite. On 
November 16, range-azimuth information was acquired by 
locating the theodolite over station MUSSEL and initialing 
cn USE MCN. The initial direction was checked by sighting 
on KM3Y MAST. Azimuth-azimuth positions were acquired cn 
November 23. A theodclite was set over USE MON and an 
initial direction was to MUSSEI. A second theodolite was 
set at MUSSEL using USE MON for the initial direction. 
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Range information was recorded using a Racal Decca 
Trispcnder system, a microwave system commonly used for 
nearshore, line-of- sight survey work. On October 28 and 
November 30, range-range data were recorded by setting 
remote units over stations BEACH LAB and MUSSEL. Before and 
after the survey, the ranging system was calibrated ever the 
fixed base line USE MCN to MUSSEL. Daily checks in the 
survey area were made to determine if the system was working 
properly. This was accomplished by maneuvering the survey 
vessel tc a point where two known navigational ranges inter- 
sected. One navigational range was formed by stations 
MONTEREI AMERICAN CAN COMPANY STACK and MONTEREY RADIO 
STATION KMBY MAST. A second navigational range was formed 
by stations MONTEREY HARBOR LIGHT 6 and MONTEREY BLUE 
LIGHTHOUSE. 

Track control for range-azimuth and range-range posi- 
tions was accomplished by steering the vessel along range 
arcs. The spacing between range arcs for most lines was 
planned to be 40 meters. Distance between positions along a 
sounding line averaged approximately 200 meters. The 
azimuth-azimuth lines were controlled by steering a magnetic 
compass heading. 

The data acquired under training conditions contained 
several deficiencies that would normally not be tolerated. 
For example, the quality of the line steering was generally 
poor; the vessel wandered off the arc more than 10 meters in 
several instances. The quality of the sounding lines run 
using azimuth-azimuth control was extremely deficient; the 
position plot of these lines show a jagged path by the 
vessel. Under normal hydrographic procedures, these posi- 
tions would be rejected. Since the intent of this study is 
to demonstrate accuracy analysis techniques, these deficien- 
cies prove to be inconsequential; the acquired data are 
adequate to demonstrate the concepts. 
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IV, BI50LTS AM data ANA LYSIS 



A. EATA PEOCESSING 

Autcmated processing of the positional survey data was 
done on the NPS IBM 370/3033AP computer system. Graphic 
displays were constructed using the Display Integrated 
Software System and Plotting Language (DISSPLA) developed by 
the Integrated Software Systems Corporation (ISSCO) 

£Bef. 16]. All computer programs involved in data 
processing were written in the WATFIV programming language. 

Ccmputations were made in an X-Y coordinate system tased 
on a Modified Transverse Mercator (MTM) projection. A MTK 
projection is essentially the same as a Universal Transverse 
Mercator (DTM) projection, the only difference being that in 
a MTM projection a central meridian is picked near the 
survey area instead of being fixed at a particular meridian 
[fief. 17], 

The central meridian, controlling latitude, and false 
easting values define the coordinate system used for compu- 
tations. The central meridian for the projection was chosen 
to be longitude 121° 52* 30" W which is approximately the 
mean longitude of the survey area. The controlling lati- 
tude, the distance in meters from the equator to a reference 
latitude, was chosen to be 4,050,000 meters. A false 
easting cf 5,000 meters was chosen as the value of the 
X-coordinate at the central meridian. 

Three shore control stations were used in the acquisi- 
tion of survey data. The geodetic positions of these 
stations were converted to the X-Y coordinate system (Table 
V) using program OCOMPS, which is a hydrographic utility 
package available to students at NPS. 
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TABLE ? 

Coordinates of Control Stations 



RATION NAME 


GEODETIC 


CCORD. 






MTM COOED 


• 


DSE MON 


360 


36' 


04.685" 


N 


Y 


= 1982.43 


m. 




1210 


52' 


35.900" 


W 


X 


= 4853.36 


IT . 


KDSSEL 


360 


37' 


18.151" 


N 


Y 


= 4247.42 


ID« 




12 10 


54' 


11.628" 


w 


X 


= 2474.75 


ID • 


BEACH LAB 


360 


36' 


05.571" 




Y 


= 2009. 86 


m. 




1210 


52' 


33.427" 




X 


= 4914.75 


m . 



B. ACCDEACY ANALYSIS OF HYDfiOGBAPHIC POSITIONING DATA 

The objective of this section is to illustrate how the 
accuracy of hydrographic positioning data can be classified 
using Burt's method of circles of equivalent probability. 
The radius of the 90 percent confidence circle was computed 
for each position; it provides a quantitative measure of 
repeatable accuracy. 

For subsequent accuracy computations, the following 
assumptions were made: 

i. The standard error for the microwave ranging 
system used in the range-range and range- 
azimuth computations is 3 meters, 

ii . For azimuth-azimuth and range-azimuth positions, 
the pointing error of the theodolite is 1.3 
meters at all ranges. 

iii. The two LOP's involved in all types 

of positioning are independent { p = 0) . 

1 2 

iv. The data are free of systematic errors. 

Eaw range and azimuth data were hand logged into a data 
file for processing. A modification of program UCOMPS was 
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used to compute X-Y coordinates of all positions. Based on 
geometric relationships discussed earlier, angles of inter- 
section of the LCP’s vere then computed for range-range and 
azimuth-azimuth points. The angles of intersection for all 
range-azimuth positicrs are 90°. 

The range-range and azimuth-azimuth data were then 
passed to WATFIV subroutine PECB (Appendix A). As input 
parameters, the subroutine accepts two standard errors of 
the TCP's and the corresponding angle of intersection. The 
output parameters include the semi-major and semi-minor axes 
of the 90 percent confidence ellipse, the radius of the 90 
percent confidence circle, and the areas covered by both 
figures. 

Subroutine PEOB uses a linear approximation to determine 
the value of the function K for varying values of the eccen- 
tricity, c, in Burt's method. A linear interpolation was 
performed by first taking the eleven discrete values of c 
and K for a probability of 90 percent from Table IV and then 
constructing a series of relationships for K as a function 
of c (Table VI) . 

Values of the radii of 90 percent confidence circles for 
range-range data were plotted at their respective positions 
(Fig. 4.1). The arcs of circles connecting the two control 
staticns BEACH LAB and MUSSEL represent lines of constant 
intersection angle (30°). Of the range-range data set, 
position 848 (Appendix B) — coordinates X = 4119.01, Y = 

4735. C7 — was found to have the smallest radius (strongest 
position) of 6.4 meters and an angle of intersection cf 
90.2°. Position 137 — coordinates X = 3345.86, Y = 

3873-34 — represents the weakest position with radius value 
of 15.3 neters and an angle of intersection of 26. 7°. 

The positional accuracy degrades rapidly as the inter- 
section angle approaches 30°; the 30° arc represents a line 
of constant 13.7 meter radius. Within 400 meters of the 30° 
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TABLE VI 

Linear ApproxiBations for K as a Function of c 






In ter val of c Linear Inter pol at ion 

~ funct ion fo r K 



0.0 


— 


0. 1 


K 




.0306c 


+ 


1.64485 


0.1 


— 


0.2 


K 


= 


-0940c 


+ 


1.63851 


0.2 


— 


0. 3 


K 


= 


. 1652c 


+ 


1-62427 


0.3 


— 


0.4 


K 


= 


.2535c 




1.59773 


0.4 


— 


0.5 


K 




.3790c 




1.54758 


0.5 


— 


0.6 


K 


= 


.5444c 


+ 


1.46433 


0.6 


— 


0.7 


K 


— 


.7101c 


+ 


1-36546 


0.7 


— 


0.8 


K 


= 


.8508c 


+ 


1-26697 


0.8 


— 


0.9 


K 


= 


. 9475c 




1 . 1896 1 


0.9 


- 


1.0 


K 


= 1 


.0361c 


+ 


1 .10987 



intersection arc, the radius varies between 8 and 15 meters. 
The radii values charge slowly in the vicinity of the 
minimuir value of 6.4 meters which corresponds to an angle of 
intersection of 90°. 

The radii of 90 percent confidence circles associated 
with the azimuth-azimuth positions acquired using control 
stations USE MON and P.USSEL were also plotted- at their 
respective positions (Fig. 4.2). The standard errors of the 
LOP’S are assumed to te 1.3 meters; the resulting improved 
accuracy is evident. The maximum value of the 90 percent 
confidence circle radii is 8.7 meters at position 
637 — coordinates X = 4327.25, Y = 28 18. 39--which corresponds 
to an angle of intersection of 159.8° (or in terms of the 
supplement, 20.2°). Position 682 — coordinates X = 4611.20, 

Y = 4421.29 — represents the strongest position recorded 
during the survey with a 90 percent confidence circle radius 
of 2.8 meters and an angle of intersection of 91.0°. 

Again, the rapid degradation of accuracy is noted 
approaching 8 = 150°. The arc of the 150° intersecticr, 
angle represents a corstant radius of 5-9 meters. Discrete 
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Figure 4. 1 Eange-Hange Accuracy Aaalysis 
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Figure 4.2 Azimuth-Azimuth Accuracy Analysis 
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values along the arc corfirm this qualitatively. A large 
area of strong positional accuracy surrounds the area where 
6 = 90°. Nuaerous values of 2.8 meters are present near 

the top cf the plot. 

Using the assumptions stated at the beginning of this 
section, the values for all radii of 90 percent confidence 
circles for range-azimuth positions are 5.1 meters. This 
computation was carried out in Example 2 of Chapter II. 

Since this case is trivial, the data are not displayed 
graphically. 

Positioning data were also classified based on the 
parameters of the 90 percent confidence ellipse. WATFIV 
program ELLIP (Appendix C) was used to generate the parame- 
ters cf the 90 percent confidence ellipse for range-range, 
azimuth-azimuth, and range-azimuth positioning data. The 
program was initialized by entering the coordinates of the 
control stations and standard errors of the LOP’s. The fix 
number, hydrographic position coordinates, and angle of 
intersection were then read in from a data file. Subroutine 

PROS was called to compute values for Ka and Ka . 

X y 

The angle of orientation of the major axis of the 
ellipse, measured clockwise from north, was then computed. 
For range-range and azimuth-azimuth positions, the LCP 
generated from the left control station was used as the base 
LOP. For range-azimuth positions, the LOP formed by the 
theodolite was used as the base LOP. First, the orientation 
of the base LOP in the coordinate system was determined. 

The orientation of the major axis of the error ellipse rela- 
tive to the base LOP (0) was then computed using Equation 
2.29. By adding or subtracting 9 to the orientation of the 
base LOP, the orientation of the major axis of the error 
ellipse in the coordinate system was determined. This angle 
takes on values from 0° to 180°. Appendix D consists cf the 
confidence ellipse classification scheme for range-range. 
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azimuth-azimuth aud range-azimuth data. Forty positions for 
each positioning geometry are listed for comparison to the 
classif ication scheme presented in Appendix 3. 

Appendix B lists the data by position number, X-'i coor- 
dinate, angle of intersection, and radius of the 90 percent 
confidence circle. Appendix D lists the data by position 

number, X-Y coordinate, angle of intersection, Ka # Ka r and 

X y 

angle of orientation for the 90 percent confidence ellipse. 
Ihese appendices are similar to hydrographic survey data 
bases and demonstrate accuracy classification schemes based 
on the tvo criteria. 

C. ACCDBACY PREEICTICHS 

The overall positional accuracy of a survey can be 
controlled by computing accuracy values before data acquisi- 
tion is begun. For example, if the hydrographer is using 
radii of 90 percent confidence circles as an accuracy 
criterion, the minimum allowable angle of intersection for 
two LCP’s can be computed for meeting specifications. The 
nature of the survey area may allow the flexibility to 
change system geometry to maximize accuracy dt a specific 
location or to maximize the area covered with a given accu- 
racy. By making accuracy computations before acquiring 
data, the hydrographer may also have the option of deciding 
what type of positioting system is to be used to meet accu- 
racy requirements. 

The construction of reliability contours is one method 

to display the expected positional accuracy. Reliability 

contours, lines of constant repeatable accuracy which are 

functions of the system geometry and standard errors of the 

positicnirg equipment, can be constructed about shore 

stations using the radii of 90 percent confidence circles 

criterion or the less desirable d^^ value. 

rras 
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Consider the equations that have been developed in 
Chapter II for the determination of radii of 90 percent 
confidence circles using Burt’s method. For uncorrelated 
lOP’s in a range-range or azimuth-azimuth system, the 
repeatable accuracy cf a hydrographic position is a function 
only cf the angle of intersection, assuming the standard 
errors of the LOP's are constant throughout the survey area. 
The Iccus of points which define a constant angle of inter- 
section for two LOP'S in a range-range or azimuth-azimuth 
system is a circle which passes through both control 
stations. Given the coordinates of the two control 
stations, the equations of these circles can be determined. 

Construction of reliability contours involves several 
simple trigonometric relationships (Fig. 4.3). let IB be 
the line connecting the two shore control stations L and E 
in a range-range system. The length of line LB is b. Ihe 
circle through both stations defines a line of constant 
intersection angle fcr two LOP's. The radius of the circle 
is r. The distance e is measured along the perpendicular 
bisector of the lice IE to the center of the circle at 
point 0(h,k) and is given by 

b 

® ■ 2tan6 (4.1) 

Knowing e and the radius r, the coordinates of point 0 can 
be ccmputed. The, equation of the circle is then 

* (x - h)* + (y - k)^ (4^2) 

These two equations were used to generate reliability 
contours for display on a computer graphics terminal. Using 
Burt's method, the angles of intersection of two LOP's were 
computed for discrete values of radii of 90 percent confi- 
dence circles. Reliability contours about stations BEACH 
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Figure 4.3 Ccistructiou of a Reliability Curve 

LAB and MOSSEL for a range-range system (a =a = 3 meters) 

1 2 

were constructed (Fig. 4.4). Osing Equation 4.2, X-Y coor- 
dinates were generated for points laying on different reli- 
ability circles. A curve-fitting subroutine in the CISSPIA 
library was used to generate the circles through the 
computed points. The 13-meter accuracy contour corresponds 
to an angle of intersection of 31.6°, while the T-meter 
accuracy contour corresponds to an angle of intersection of 
67.9°. The best achievable accuracy of the system is 6.4 
meters at 9Qo. 

For comparison purposes, reliability contours were 

constructed about BEACH LAB and MUSSEL for azimuth-azimuth 

geometry { o =a =1.3 meters). The increased accuracy of 
> 2 

this configuration is evident (Fig. 4.5) . 



6 1 



The 3-meter 



contour corresponds tc an angle of intersection of 69. 4c 
while the 6-ineter contour corresponds to an angle of inter- 
section of 29.6°. The best achievable accuracy at an inter- 
section angle of 90° is 2.8 meters- 

A second scheme was used tc display accuracy predictions 
for the twc positioning methods. Given the coordinates of 
ESACH lAB and MUSSEL, a series of discrete points spaced 800 
meters apart, were generated throughout the survey area. 

The values for the radii of 90 percent confidence circles 
were then computed at each point with the use of subroutine 
EROB. Figures 4.6 and 4.7 illustrate this prediction 
scheme. These figures present the same information as 
Figures 4.4 and 4.5 in a different manner. The 30° angle of 
intersection contour is shown on both figures. 
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Figure 4.4 Heliaiiility Contours: Range-Range Geometry 
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Figure 4.5 Beliability Contours: izimuth-Azimath Geometry 
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ligure 4*6 Hauge-fiange Point Accuracy Prediction 
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Figure 4-7 Azimu th-Aziauth Point Accuracy Prediction 
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V. C0N CI DS10N5 R ECOMaEM DA TIONS 

A. ACCUEACI SPECIFICATIONS 

iDterpretation of the 1982 IHO positioning standards in 
terms cf 90 percent ccnfidence circles yields some inter- 
esting results with respect to present day survey practices. 
For example, for a 1 : 10,000-scale hydrographic survey, NOS 
usually uses microwave positioning systems in a range-range 
mode, and assumes a standard error of 3 meters for each LCP. 
Surveys are frequently conducted between the 30° to 150° 
angle of intersection limits. Osing the 90 percent ccnfi- 
dence circle critericr, the radius of the circle should net 
exceed 10 meters. Hewever, the radius value for 6 = 30° and 
150° is 13.7 meters. The values of Ka and Ka for the 90 

X y 

percent confidence ellipse are 17.6 and 4.7 meters, respec- 
tively. To meet the SO percent criterion for a 1:10, 000- 
scale survey, the 6 limits should be 42° to 138°. 

A ximuth-azimuth positioning is accurate enough for 
1: 5 , 0 0 C-scale surveys, using 6 limits of 35° to 145°, 
assuming a standard error of 1.3 meters for each LCP. liith 
the standard error assumptions used for range-azimuth, the 
90 percent radius is 5.1 meters for all positions. Given 
the uncertainties of the standard error figures, it is 
rational to assume that range-azimuth positions can meet the 
5-meter accuracy standard for 1: 5, 000-scale surveys. In 
fact, range-azimuth positional accuracy can exceed azimuth- 
azimuth accuracy when the later's 6 is less than 35°. For a 
3-meter a range-range configuration, it is impossible to 
meet 1:5,000 specifications with any 6- 

As a general guideline, the 30° to 150° angle of inter- 
section limit is a good rule to use for uncorrelated LOF's. 
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However, as mentioned for 1 : 10 ,000-scale surveys in a range- 
range mode {0= 3 meters), this rule does not always hold. 

On the other hand, it is possible to have g 's of less than 
30° and still meet specif icaticns. For example, azimuth- 
azimuth positioning can theoretically be used for 6’s of 18° 
to 162° for a 1: 10,OOC-scale survey. However, the eccen- 
tricity cf the error ellipse is so small that the distortion 
introduced by using confidence circles can become 
misleading. In view cf this, eccentrici ties of less that 
0.2 should not be used. 

Using the 90 percent radius criterion, a table has been 
assembled illustrating the 6 limit for various positioning 
geometries at different survey scales, using assumed stan- 
dard errors (Table VII) . The information in Table VII 
illustrates that the 30° to 150° 6 limit need not be fixed. 
The 6 limits should vary based on the scale of the survey 
and the precision of the positioning equipment. Accuracy 







TABLE 7II 






g Limits for Surveys 




Survey 

Scale 


90% 

Radius 




R-R 

io = 10) 
B Limit 


Az- Az 
(0 = 1.3) 
6 Limit 


1:2,500 
1:5,000 
1 : 10,000 
1:20,000 i 
1 :h0,000 


(m) 

2.5 I 

5. 0 
10.0 
20.0 
40. 0 


(deg) 


(deg) 


(deg) 






35-145 

23-157* 

23-157* 

23-157* 


42-138 
27- 153 
23-157* ] 


35-145 


* Eccentricity limit of 0.2 
Note: 90% radii cf all range- 

assumed to be 5. 1 meters for 


-azimuth positions are 
0=3 and a = 1.3. 
1 2 
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figures as a function of B for uncorrelated LCP’s have teen 
compiled using standard errors of 1.3 meters for azimuth- 
azimuth (Table VIII) ,3.0 meters for range-range short-range 
(Table IX), and 10 meters for range-range medium-range 
(Table X) positioning systems. 



TABLE 7III 

Accuracy Figures for o =o =1.3 m, o =0 

12 12 



Angle of Ko 


K 0 


Radius of 


Area of 


Area of 


Infer. 


X 


y 


905 Circ. 


Ellipse 


Circle 


(deg) 


(m) 


(m) 


(m) 


(sg m) 


(sg m) 


90 


2.8 1 


2.8 


2. 8 j 


24 


24 


85 


2.9 


2.7 


2.8 


25 


25 


80 


3. 1 


2.6 


2.8 


25 


25 


75 


3.2 


2. 5 


2.9 


25 


26 


70 


3.4 


2.4 1 


3.0 


26 


28 


65 


3.7 


2.4 


3.1 


27 


30 


60 


3. 9 


2. 3 


3.3 


28 


34 


55 


4.7 


2.3 


3.5 


30 


38 


50 


5. 1 


2. 2 


3.8 ] 


32 


44 


45 . 


5. 1 


2. 1 


4. 1 


35 


53 


40 


5.8 


2. 1 


4.5 


38 


65 


35 


6.6 


2. 1 


5. 1 


43 


83 


30 


7.6 


2.0 


5.9 


49 


1 10 


25 


9. 1 


2.0 


7.1 


58 


156 


20 


11.4 


2.0 


8. 8 


71 


24 1 


15 


15. 1 


2.0 


11.6 


94 


425 


10 1 


22.6 


2.0 1 


17.4 1 


14 1 1 


949 


5 


45. 2 


2.0 1 


34.7 


281 


3, 78 1 



K = 2.146 for 905 probability 
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TABLE X 

Accuracy Figures for .a = a = 10 p = 0 

12 12 



Angle of Ka 


Ka 


Radius of 


Area of 


Area of 


Inter. 


X 


y 


90% Circ. 


Ellipse 


Circle 


(deg) 


(m) 


(n) 


(n) 


(sg m) 


(sg ffl) 


90 


21.5 


21.5 


21.5 


1,447 


1,447 


85 


22.5 


20.6 


21.6 


1, 452 


1,45 5 


80 


23.6 


19.8 


, 21.8 


1,469 


1,497 


75 


24.9 


19. 1 


22.3 


1, 49 8 


1,562 


70 


24.5 


18.5 


23.0 


1, 540 


1,657 


65 


28.2 


18.0 


23.9 


1,596 


1,79 8 


60 


30.3 


17.5 


25.2 


1,671 


1,989 


55 


32.9 


17. 1 


26.8 


1 , 766 


2,252 


50 


35.9 


16.7 


28.9 


1,889 


2,615 


45 


39.6 


16.4 


31.5 


2, 046 


3,117 


40 


44.4 


16.1 


34.9 


2^251 


3"835 


35 


50.5 


15.9 


39.5 


2, 522 


4,890 


30 


58.6 


15.7 


45.6 


2,894 


6,528 


25 


70.1 


15.5 


54.3 


3, 423 


9,250 


20 


87.4 


15.4 


67.4 


4,230 


14,270 


15 


116.3 


15.3 


89.4 


5,590 


25, 126 


10 


174. 1 


15.2 


133.7 


8,332 


56^130 


5 


347. 9 


15.2 


266.9 


16, 600 


223,727 



K = 2.146 for 901 probability 
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B. OSES FOE ACCOEACl FIGORES 



NCS is currently developing the Shipboard Data System 
III (SDS III) , a hydrographic data acquisition and 
processing system which will replace the present HYDROLCG/ 
HYDRCELOI system. SDS III will revolutionize data acquisi- 
tion and processing techniques with the capability to 
perfori high-speed calculations and display color graphics, 
^iith this increa-sed computer potential, data manipulat icns-- 
such as accuracy computations--can be performed. 

Each position in a survey can be given a quality figure 
based on the radius of the 90 percent confidence circle. 

This figure is sufficient for non-critical positions of 
ordinary hydrographic data. Critical positions are those 
which are determined for significant features (i.e., wrecks, 
least depths, rocks, and other potential hazards). For 
these positions, the parameters of the 90 percent error 
ellipse can be computed, as well as the radius of the 90 
percent' confidence circle. 

Many schemes can be envisioned for the use of an accu- 
racy figure. For example, suppose the position of a 
submerged pile was determined by range-azimuth geometry in a 
prior survey. The radius of the 90 percent confidence 
circle is then 5.1 meters (Ex. 2, Ch. II). The charting 
agency now wishes to relocate the pile to determine if it 
still exists and is still a hazard to navigation. In low 
water visibility, a common technique used to resolve such an 
item would be to send divers down over the reported position 
and conduct a circle search. One diver remains at the 
reported position, holding a line, while the other diver 
swims a circumference holding the other end of the line. 
Theoretically, if the line is about 5 meters long and a hang 
does not occur, it is 90 percent certain that the pile has 
been removed. For a higher confidence, the line is 



lengthened. In an investigation such as this, it is advi- 
sable tc be conservative and use the maxiinuin length of line 
which is oferationall j feasible to provide coverage cf an 
area as large as possible. The radius of the 90 percent 
confidence circle gives the hydrographer a rough figure fcr 
answering the guesticn: Does the submerged pile exist? 

Knowing the paraieters of the error ellipse could be 
useful fcr conducting wire-drag, wire-sweep, and side scan 
sonar operations. Fci a position obtained with low preci- 
sion positioning equipment, the search to relocate a 
submerged feature could cover a large area. Knowing the 
parameters cf the error ellipse could reduce the area, time, 
and effort of the search. The search pattern could be 
planned to cover the desired confidence ellipse. 

With the quantification of accuracy, a decision must he 
made concerning how much confidence is needed to delete a 
certain feature from the chart after a search has teen made. 
The 90 percent confidence level may be too low, whereas tte 
S5 or 99 percent level may suffice. A balance must te main- 
tained between confidence of disproval and time and effort 
spent on the search. 

Accuracy predictions in the form of reliability contours 
can be displayed using computer graphic terminals. These 
displays will contribute to the efficient planning of 
surveys to meet specif ications. Given the survey area, the 
available control, the positioning methods, and the preci- 
sion of the positioning equipment, the hydrographer can plan 
the accuracy of the survey before it is conducted. The 
survey area and the available control may be such that there 
is flexibility to change control stations to optimize accu- 
racy over an area of critical importance. This information 
can be displayed graphically and plans for the survey can be 
made accordingly. Likewise, given an accuracy limit, such 
as a 10-meter radius cf the 90 percent confidence circle, 
the area tc be covered at that accuracy can be maximized. 
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Many variables eiist when considering accuracy reguire- 
ments fcr a hydrographic survey. In general, higher accu- 
racy means more time, money, and effort. Azimuth-a 2 imuth 
geometry is the most accurate method of positioning analyzed 
in this thesis. This method involves at least two people 
ashore and good ship-to-shore communications. Currently, 

NOS acquires these data manually, which minimizes the speed 
that the vessel can operate and adds to processing time. Or. 
the other hand, a survey using a medium-range system needs 
little shore support and the data acquisition is automated. 
Accuracy predictions help keep a balance between accuracy 
and effort. If the desired accuracy is attainable using a 
range-range system instead of an azimuth-azimuth system, 
then the choice is obvious. 

Hydrographic positioning in the future will be dominated 
by two methods. For offshore surveys, the Global 
Positioning System (GPS) is expected to give positional 
accuracy to 10 meters or better. GPS is a satellite posi- 
tioning system currently being deployed by the Department of 
Defense and will provide near worldwide coverage for users. 
Since the full constellation of 18 satellites will not' be 
operational until 1988, it is not yet known if the expected 
accuracy of 10 meters will be met. Nearshore surveys may 
use multiple LOP’S for establishing hydrographic positions. 
The principle of least squares is applied to redundant 
observations yielding the most probable position. Fcr both 
GPS and least squares positioning, confidence ellipses and 
circles can be determined, although the techniques involved 
are much mere complicated than those presented in this 
thesis. 

The accuracy classif ication scheme presented in this 
thesis is predicated on the elimination of systematic 
errors. Much work is needed in identifying the sources of 
systematic errors associated with hydrographic positioning 
equipment. 
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APPENDIX A 

SOEBCOTINE FOB 90 FEBCENT CONFIDENCE CIRCLE PARAMETEES 



SUEECUTINE PEOE (SI G 1 , S IG2 ,C0R , TBET A, SGX90, SGY90, 

* RADIOS, ELAR.CIBAE) 

IMIIICIT EEAL*i» (A-H-O-Z) 

COMPOTES RADIOS OF 90% CONFIDENCE CIRCLE (BCET, METHOD 2) 
THIS SOBEOUTINE WORKS FOR CORRELATED AND ONCOERELATED 
LINES OF POSITION. 

INEDT lABAMETEES: 

SIG1 AND SIG2- STANDARD ERRORS OF TWO LOP’S 
IBETA - ANGLE OF INTERSECTION IN DEGREES 

(0-180 DEGl 

COB - CORRELATION COEFFICIENT 

(USUALLY ZERO EXCEPT FOR HYPEEECIIC 
OR SEXTANT POSITIONING) 



OUTPUT PARAMETERS: 

SGX90 AND SGYSO- 

EADIUS 

ELAR 

CIRAR 



SEMI-MAJOR AND MINOR AXES 
OF 90^ ERROR ELLIPSE 
RADIUS OF 90% CONFIDENCE CIRCLE 
AREA OF 90% CONFIDENCE ELLIPSE 
AREA OF 90% CONFIDENCE CIRCLE 



WORK KITH AN ANGLE lESS THAN 90 
IF (TEETA .GT.90.) BETA =180. 

IF TBETA.LE.90.) BETA=TBETA 
CHANGE DEGREES TO RADIANS 
RAD=.0174532*BETA 
TRAKSFCRHATION SIGl AND SIG2 TO 
RECTANGULAR SYSTEM 

SIGA=SQRT ( (1./ ( (SIN (RAD) ) = 

=»SIG2*COS (RAD) +SIG2**2) -SIC 
SIGB=SIG2 

TRANSFORM CORREIATICN COEFFICIENT 
RECTANGULAR COORDINATE SYSTEM 

A=i(SIG2*COS (EAE) )/SIG1) +COR 

F = 1/SQRT ( 1+2*CCE*SIG2*COS (RAD) /SIG 1+ (SIG2/SIG1) **2* 
* (CCS (EAD)^) **2) 

COEAE=A*F 



DEGREES 

,-TBETA 



CORRELATED 



** 2 ) )* (SIG1**2+ 

• 2 ] 



CG2**2) 



2. *COR*SIG 1* 



TO CORRELATED 



TRANSFORM TO UNCORRELATED RECTANGULAR 
AA=SQRT { fSIGA**2+SIGB**2)/2) 

CC=SQRT (1- (4*S1GA**2*SIGB**2* (1-CORAB**2) )/ 
* (SIGA=»*2 + SIGB**2) **2) 

DD=SQRT (1+CC) 

SIGX= AA*DD 

SIGY=SQRT (SIGA**2+SIGB**2-SIGX**2) 



COMPUTE ECCENTRICITY OF ELLIPSE 
C=SIGY/SIGX 

C COMPUTE BURT’S K FACTOR BY LINEAR INTERPOLATION 
IF J[ (C. LE. 1 . ) .AND. jfC.GT. 0.9) ) THEN 
B=T.0361*C+ 1.10987 

ELSE IF ((C.LE . 0 . 9) . A ND. ( C. GT. 0 . 8) ) THEN 
B=0.9475*C+1. 1896 1 

ELSE IF ((C.IE.0.8) .AND. (C.GT.0.7)) THEN 
B=0.8508*C + 1. 26697 

ELSE IF ( (C.IE. 0.7) .AND. (C. GT. 0. 6) ) THEN 
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B=0.7 10 1*C + 1.36546 

ELSE IF { (C .IE. 0.6) .AND. (C.GT. 0.5) ) THEN 
B=C. 5444*0+1.46488 

EISE IF ({C. IE. 0.5) .AND. (C.GT. 0.4)) THEN 
8=0.3790*6+1.54758 

ELSE IF ( (C. IE. 0.4) .AND. (C.GT. 0.3) ) THEN 
8=0.2535*6+ 1.59778 

ELSE IF (fC. IE. 0.3) .AND. (C.GT. 0.2)) THEN 
8=0.1652*6+1.62427 

EISE IF (Ic. IE. 0.2) .AND. (C.GT. 0. 1)) THEN 
B=.094*C+1. 63851 



ELSE 

8=0. 0306*C+1. 64485 

END IF 

EADIDS=B*SIGX 
SGX90 = 2. 146*SIGX 
SGY90=2.146*SIGY 
CIEAR=3. 1415S26*EADIUS**2 
E1AE=3. 1415926*SGX90*SGY90 
EETOEN 
END 
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APPENDIX B 



ACCOBACI CLASSIFICAIION: 90 PERCENT CONFIDENCE CIRCLES 

CLASSIFIED RANGE-RANGE POSITIONS 
Control Stations: BEACH LAB 1982 and MUSSEL 1932 

Standard Error Used in Computations: 3 meters 



‘ix 


X 


Y 


Angle of 


Radius 


[2- 


Coordinate 


Coordinate 


Intersection 


90^ Ci. 


1 


2668. 05 


4942.07 


127. 0 


8.2 


2 


2852. 01 


5076. 69 


121. 6 


7.7 


3 


3041. 58 


5194.85 


113. 6 


7.4 


4 


3040. 27 


5148. 74 


1 17. 0 


7.3 


5 


2838. 26 


5021.35 


120.2 


7.6 


6 


2640. 06 


4877. 13 


126.9 


8.2 


7 


2553.46 


4752.36 


130.4 


8. 7 


8 


2724. 10 


4885.47 


121.3 


7.7 


9 


2889.06 


5013.63 


1 17. 6 


7.3 


1 0 


3075.61 


5124. 37 


115.0 


7.2 


1 1 


3172. 63 


5136. 23 


1 12.7 


7.0 


12 


2958. 13 


5006. 65 


1 14. 4 


7. 1 


13 


2771. 33 


4876. 02 


117.9 


7.4 


14 


2581. 02 


4729. 89 


126. 9 


8.2 


15 


2584.4 1 


4665.31 


124. 0 


7.9 


16 


2740. 05 


4805.63 


116.7 


7.3 


17 


2913. 27 


4928.51 


112.8 


7.0 


18 


3097.44 


5047. 77 


111.2 


6.9 


19 


3193. 94 


5103. 46 


110.9 


6.9 


20 


2627.47 


4659. 05 


118.8 


7. 4 


2 1 


2762. 15 


4769.60 


113.2 


7. 1 


22 


2904. 91 


4875.24 


110.5 


6.9 


23 


3056.38 


4975.73 


109.3 


6.9 


24 


3207.36 


5064. 00 


108. 9 


6. 8 


25 


3373.36 


5101. 96 


107. 1 


6.8 


26 


3190. 64 


5007. 38 


106. 8 


6.8 


27 


3015. 45 


4900. 75 


107. 1 


6.8 


28 


2839.38 


4778. 31 


108. 7 


6.8 


29 


2679.99 


4651.31 


112.8 


7.0 


30 


2727.88 


4632. 77 


106. 9 


6.8 


3 1 


2899.83 


4776.97 


105.2 


6.7 


32 


3092.77 


4901. 75 


104.4 


6. 6 


33 


3295.42 


5019. 92 


105.0 


6.7 


34 


3502. 07 


5123. 36 


106. 0 


6.7 


35 


3697.02 


5151.09 


105.3 


6.7 


36 


3474. 77 


5065. 57 


104.0 


6. 7 


37 


3257.48 


4953. 37 


102.7 


6.6 


38 


3043.03 


4821. 88 


101.7 


6.6 


39 


2845- 07 


4680.80 


101.7 


6.6 


40 


2746. 96 


4608. 32 


103. 1 


6.6 


41 


2748. 09 


4550. 40 


97.5 


6.5 


42 


2931.78 


4701.21 


98.4 


6.5 


43 


3134. 33 


4834.56 


99.4 


6. 5 


44 


3317. 94 


4944.55 


101.0 


6.6 


45 


3515. 14 


5043. 56 


102.6 


6.6 


46 


3720. 18 


5117. 15 


103.9 


6.6 


48 


3724. 93 


5083. 62 


102.6 


6.6 


49 


3514.33 


4989. 13 


100.3 


6.6 


50 


3312. 88 


4891. 55 


98.5 


6.5 
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Fix 

No. 

51 

52 

53 

54 

55 

56 

57 

58 

59 

6 1 

62 

62 

64 

65 

66 

67 

68 

69 

70 

7 1 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

84 

84 

85 

86 

87 

88 

89 

90 

9 1 

92 

93 

94 

95 

96 

97 

98 

100 

101 

102 

103 

104 

105 

106 

107 

108 

110 

1 1 1 

112 

113 

114 

115 

116 

1 17 



RANGE-RAKGE ACCURACIES (CONTINUED) 



X 


Y 


Angle of 


Radius of 


Coordinate 


Coordinate 


Intersection 


90^ Circle 


3114.26 


4776.58 


96.5 


6.5 


2921. 85 


4647. 63 


94.8 


6. 5 


2753.76 


4502. 52 


91 . 5 


6.4 


2793.3 1 


4487.76 


86.6 


6. 4 


2960. 81 


4619.38 


90.6 


6. 4 


3160. 14 


4750. 25 


93.6 


6.4 


3325. 1 1 


4847.99 


96.0 


6.5 


3492. 32 


4942.05 


98.4 


6. 5 


3693. 83 


5032. 60 


100.8 


6.6 


3568. 42 


4910.65 


96.3 


6.5 


3370.43 


4812. 54 


93.4 


6.4 


3174.72 


4695. 4 1 


89.7 


6.4 


2998.35 


4566.29 


84.5 


6. 5 


2824.36 


4431, 37 


76.9 


6.6 


2819. 19 


4368.35 


67. 7 


7.0 


3002. 82 


4519. 53 


80.0 


6.5 


3184. 74 


4657. 79 


86. 9 


6. 4 


3378. 02 


4763.00 


90.5 


6. 4 


3578. 71 


4871.32 


94. 4 


6.5 


3788. 85 


4962. 48 


97.7 


6.5 


4010. 25 


5031.39 


100. 4 


6.6 


3732. 84 


4895. 17 


95.0 


6.5 


3509.99 


4799.08 


91.3 


6.4 


3307. 18 


4688. 48 


86.9 


6, 4 


3094.50 


4546. 29 


80. 1 


6.5 


2901. 06 


439 1. 68 


68.5 


7.0 


2754.33 


4263. 59 


49. 5 


8.7 


3754. 15 


5020.93 


100. 1 


6.5 


3514.53 


4919. 46 


97. 2 


6.5 


3287.25 


4788.27 


93.3 


6.4 


3062.52 


4651. 12 


89.4 


6.4 


2846. 46 


4482. 96 


82.4 


6.5 


2845. 13 


4562. 17 


91.3 


6.4 


2995. 88 


4403, 69 


68.0 


7.0 


3165. 76 


4254,99 


52. 7 


8.3 


3331. 61 


4102.46 


43.3 


9.8 


3500.54 


3946. 44 


37. 5 


11.1 


3672. 90 


3785.70 


34. 0 


12. 2 


3846.42 


3623. 1 9 


32.0 


12,9 


4022. 93 


3464. 18 


31.6 


13.0 


4406. 76 


3998. 00 


68. 3 


7.0 


4263.63 


4159. 42 


70.3 


6.9 


4112. 66 


4331.20 


73.9 


6.7 


3959. 57 


4527. 77 


79.9 


6.5 


3810. 03 


4681.47 


85.5 


6. 5 


3672.32 


4864. 18 


93.7 


6.5 


3535.90 


5040. 98 


102. 3 


6.6 


3065. 44 


4352. 97 


61.8 


7.4 


3068.28 


4166.32 


41.6 


10.1 


3077. 08 


4046. 24 


29.5 


13.9 


31 17. 61 


4262. 88 


52. 8 


8.3 


3062. 15 


4499. 1 1 


76.5 


6.6 


3011. 82 


4658.04 


91.7 


6.4 


3122. 9 1 


4453. 24 


71 . 4 


6.8 


3152. 10 


4199.54 


47. 1 


9. 1 


3118. 73 


4025. 84 


29.3 


14.0 


3188.79 


4154. 61 


43. 8 


9.7 


3177.20 


4396.30 


65.9 


7. 1 


3111.79 


4583. 07 


82.8 


6.5 


3182. 67 


4534. 87 


77.6 


6.6 


3099.29 


4691.03 


91.3 


6.4 


3209.60 


4484.20 


73.3 


6.7 


3235.42 


4271. 54 


55.2 


8.0 


52C4. 85 


4043. 52 


34.3 


12.0 
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RANGE-EflNGE ACCURACIES (CONTINUED) 



Fix 


X 


Y 


Angle of 


Radius of 


No, 


Coordinate 


Coordinate 


Intersection 


90S Circle 


120 


3261. 07 


411 1. 17 


42.0 


10.0 


121 


3262. 51 


4368.76 


63.7 


7.3 


122 


3175.94 


4613.63 


83.8 


6.5 


123 


3 189.26 


4667.90 


87.5 


6.4 


124 


3286.57 


4471.03 


71.9 


6.8 


125 


3317. 76 


4240. 01 


53.9 


8. 1 


126 


3282. 62 


4020. 86 


35. 3 


1 1.7 


127 


3251.48 


3951. 93 


28.6 


14.3 


128 


3309.62 


3962. 87 


31. 8 


13.0 • 


129 


3362. 57 


4200. 72 


51.7 


8.4 


130 


3331.06 


4454. 74 


70.7 


6. 9 


13 1 


3252.99 


4654.02 


85.4 


6.5 


132 


3286.69 


4673. 74 


86.3 


6.5 


133 


3366.26 


4474.82 


72.2 


6.8 


134 


3387.38 


4275.0 1 


57.7 


7. 8 


135 


3376.70 


4076.27 


42.6 


9.9 


136 


3337.33 


3910.52 


29. 0 


1 4. 1 


137 


3345. 86 


3873.34 


26.7 


15.3 


138 


3422.73 


4077.78 


44.0 


9.6 


139 


3431. 43 


4314. 80 


61.3 


7.4 


140 


3396.21 


4530. 82 


76.0 


6.6 


14 1 


3326. 96 


4700.26 


87.4 


6.4 


142 


3361.30 


4694. 84 


86.7 


6.4 


143 


3436. 64 


4506. 04 


74.4 


6.7 


1 44 


3474.89 


4309. 23 


61.5 


7.4 


145 


3464. 73 


4120.80 


48.2 


8. 9 


146 


3425.35 


3933. 68 


34. 0 


12.2 


147 


3459.34 


3928.23 


34.8 


11.9 


148 


3513.95 


4137. 04 


50. 6 


8.5 


149 


3512. 94 


4390. 84 


67.4 


7. 0 


150 


3473.67 


4567. 85 


78.4 


6.6 


151 


3385. 01 


4766. 90 


90.7 


6.4 


152 


3390.09 


4801. 65 


92. 5 


6.4 


153 


3472. 78 


4657.67 


83.8 


6.5 


154 


3527.62 


4492. 63 


73.9 


6.7 


155 


3555. 04 


4377. 14 


67.0 


7. 1 


156 


3561. 04 


4288.00 


61.4 


7.4 


157 


3543. 95 


4204.99 


55.7 


7.9 


158 


3543.69 


4262.39 


59.5 


7. 6 


159 


3543. 85 


4068. 09 


46.8 


9. 1 


160 


3493.62 


3880. 1 8 


32.9 


12.5 


719 


3519. 28 


3843.16 


31.6 


13.0 


720 


3586.31 


4061.28 


47.6 


9.0 


72 1 


3589.75 


4237. 59 


58.7 


7-7 


722 


3568.78 


4441. 1 5 


71. 1 


6.8 


723 


3525.50 


4623. 58 


81.7 


6.5 


724 


3442.67 


4796. 98 


91.7 


6.4 


725 


3478.48 


4818. 98 


92.6 


6.4 


726 


3558.35 


4658. 77 


83. 7 


6.5 


727 


3612. 48 


4473.89 


73.4 


6.7 


728 


3635.37 


4279.04 


62. 1 


7.4 


729 


3618. 98 


4069. 08 


49.0 


8.8 


730 


3574.84 


3866. 91 


35. 1 


11.8 


731 


3587. 72 


3809.47 


32. 1 


12.8 


732 


3613.95 


3910.30 


39. 1 


10.7 


733 


3648.48 


4007.29 


46.1 


9-3 


73 4 


3676.20 


4121. 4 1 


53. 6 


8.2 


735 


3677.60 


4225. 91 


59.8 


7.6 


736 


3676.79 


4332. 54 


66.0 


7. 1 


737 


3658. 87 


4436.08 


71.7 


6.8 


738 


3638. 19 


4545.66 


77. 6 


6.6 


739 


3604.38 


4649.50 


83.2 


6.5 


740 


3572.48 


4754. 1 8 


88.7 


6.4 


74 1 


3514. 26 


4836.98 


93.2 


6.4 
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RANGE-EAKGZ ACCUESCIES (CONTINUED) 



Fix 


X 


I 


Angle of 


Radius of 


No- 


Coordinate 


Coordinate 


Intersection 


S0% Circle 


742 


3570. 50 


4836. 50 


92.8 


6.4 


743 


3618. 79 


4723. 26 


87.0 


6.4 


744 


3662. 90 


4602. 44 


80.9 


6. 5 


745 


3691. 10 


4478.20 


74.4 


6.7 


746 


3715. 24 


4341.39 


67. 1 


7.0 


747 


3720. 11 


4205. 94 


59.5 


7.0 


748 


3703. 21 


4071. 17 


51.4 


8.5 


749 


3678. 03 


3939. 55 


43. 0 


9.8 


750 


3636. 04 


3801. 29 


33.7 


1 2. 3 


751 


3630. 00 


3715. 02 


28.3 


14.5 


752 


3678. 77 


3812. 88 


35.7 


11.6 


753 


3715.73 


3925.56 


43. 4 


9.7 


754 


3731. 48 


4050. 62 


51.0 


8.5 


755 


3756.38 


4168. 71 


58. 3 


7.7 


756 


3754. 68 


4300. 89 


65.5 


7. 1 


757 


3743. 02 


4427. 92 


72.2 


6.8 


758 


3710. 29 


4554. 17 


78.6 


6.6 


759 


3680.62 


4673. 88 


84.6 


6.5 


760 


3627.42 


4792.51 


90.5 


6.4 


76 1 


3572.80 


4899. 46 


95. 8 


6. 5 


762 


3627. 15 


489 1.25 


95.1 


6.5 


763 


3683.33 


4778.34 


89. 7 


6.4 


764 


3725.77 


4656.01 


83.9 


6.5 


765 


3757. 98 


4525.66 


77.5 


6. 6 


766 


3781. 56 


4403. 13 


71.5 


6.8 


767 


3798.56 


4272. 52 


64. 8 


7.2 


768 


3787. 95 


4134, 55 


57. 1 


7. 8 


769 


3772. 73 


4001. 69 


49. 4 


8.7 


770 


3738. 73 


3870. 33 


41 . 1 


10.2 


771 


3695.62 


3745. 45 


32.6 


12.7 


772 


37 12. 52 


3708.43 


31.2 


13.2 


773 


3777. 45 


3829. 70 


40. 2 


10.4 


774 


3803.43 


3957. 03 


48.0 


9.0 


775 


3818.28 


4087. 24 


55. 4 


8.0 


776 


3631. 33 


4211. 74 


62.3 


7.4 


777 


3823.70 


4338. 97 


68.8 


6. 9 


778 


3813.4 1 


4476. 12 


75.6 


6.7 


779 


3777. 13 


4614.24 


82. 1 


6.5 


780 


3739. 58 


4736.39 


87.8 


6.4 


781 


3691.63 


4868. 16 


93.9 


6. 5 


782 


3706.65 


4918. 98 


96.0 


6.5 


783 


3749.46 


4806. 82 


91. 1 


6.4 


784 


3799. 89 


4705. 41 


86.6 


6. 4 


785 


3829.30 


4581. 1 4 


81.0 


6.5 


786 


3857. 30 


4455. 00 


75. 1 


6.7 


787 


3868.49 


4327. 66 


69.0 


7. 0 


788 


3672. 94 


4198. 56 


62.5 


7.3 


789 


3862.59 


4062. 95 


55.3 


8.0 


790 


3841. 17 


3928.97 


47. 7 


9.0 


791 


3804. 62 


3797. 1 4 


39.4 


1 0.6 


792 


3759.45 


3596. 86 


27. 1 


15. 1 


793 


3808. 16 


3722.44 


35. 6 


11.6 


794 


3851. 17 


3827.55 


42.7 


9.9 


795 


3879.16 


3942. 93 


49. 6 


8.7 


796 


3896. 63 


4068. 89 


56.5 


7.9 


797 


3903. 40 


4184. 15 


62. 5 


7.3 


798 


3916.60 


4306. 31 


68.8 


6.9 


799 


3904.40 


4428. 30 


74. 5 


6.7 


800 


3879.66 


4554. 38 


80. 2 


6.5 


801 


3851. 19 


4684.13 


85.9 


6. 5 


802 


3810. 16 


4796. 95 


90.7 


6.4 


803 


3752. 88 


4908. 69 


95.5 


6.5 


804 


3706. 66 


4998.96 


99.4 


6. 5 


805 


3781.90 


4945. 89 


97.0 


6. 5 
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Fix 

No. 

806 

807 

808 

809 

8 10 

811 

812 

813 

814 

815 

816 

817 

8 18 

819 

820 

821 

822 

823 

824 

825 

826 

827 

828 

829 

830 

831 

832 

833 

834 

835 

836 

837 

838 

839 

84C 

841 

842 

843 

844 

845 

846 

847 

848 

849 

850 

851 

852 

853 

854 

855 

856 



SANGE-HANGE ACCURACIES (CONTINUED) 



X 


Y 


Angle of 


Radius of 


Coordinate 


Coordinate 


Intersec ti on 


90% circle 


3630. 54 


4843. 53 


92.3 


6.4 


3877.77 


4734. 1 8 


88.3 


6.4 


3907.61 


4618. 18 


83.4 


6.5 


3932. 77 


4495. 95 


78. 1 


6.6 


3948. 19 


4376.06 


72.8 


6.8 


3956.79 


4251. 77 


67.0 


7. 0 


3953. 81 


4125.89 


60.9 


7. 5 


3938.32 


3995. 79 


54. 1 


8. 1 


3909. 94 


3869. 78 


46.9 


9. 1 


3864. 71 


3744. 84 


38. 9 


10. 7 


3816. 98 


3619. 07 


30.6 


13. 4 


3711.62 


5131.00 


104. 5 


6.7 


3788. 16 


5014.61 


99.7 


6.5 


3851.30 


4889.79 


94.7 


6. 5 


3901. 13 


4761. 08 


89.6 


6.4 


3961.73 


4590.52 


82. 7 


6.5 


39 77. 10 


4485.23 


78.2 


6.6 


3993.32 


4118.31 


61. 5 


7.4 


3996. 85 


4193. 81 


65.2 


7.2 


3981.72 


4048. 49 


57. 9 


7.7 


3956. 88 


3905.60 


50. 2 


8.6 


3916.26 


3758.91 


41.5 


10. 1 


3826.04 


3553. 55 


27.6 


14. 8 


3961.88 


4691, 6 1 


87. 1 


6.4 


3907. 91 


4827.37 


92.4 


6.4 


3841.30 


4953.23 


97.3 


6.5 


3789. 42 


5085. 46 


102.4 


6.6 


3772.12 


5177.05 


105. 8 


6.7 


3843.76 


5069. 86 


101.7 


6.6 


3906. 65 


4961.46 


97.6 


6.5 


3959.36 


4844. 18 


93.3 


6.4 


4002. 25 


4720.39 


88.6 


6.4 


4057. 98 


4659. 81 


86.7 


6.4 


3992.26 


4828.76 


92.8 


6. 4 • 


3954. 08 


49 15. 16 


96.0 


6. 5 


3905. 75 


50 46. 6 3 


100. 8 


6.6 


3822. 81 


5167.39 


105. 2 


6.7 


3844. 43 


5129. 36 


103. 8 


6.6 


3958. 55 


5033. 27 


100.3 


6.6 


4016.73 


49 19. 44 


96. 4 


6.5 


4061.27 


4797.13 


92. 1 


6.4 


4159. 72 


4597.99 


35.5 


6.5 


41 19. 01 


4735.07 


90.2 


6.4 


4077. 07 


4872. 84 


95.0 


6.5 


4021.39 


4990. 89 


99.0 


6.5 


3961.68 


5118. 03 


103. 3 


6.6 


3992. 88 


5139.94 


104.0 


6.7 


4056.51 


5026. 15 


100.3 


6.6 


4 121. 60 


4874.08 


95.4 


6. 5 


4154. 57 


4779.44 


92.2 


6. 4 


4196.42 


4617. 75 


86.7 


6.4 
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CLASSIFIED AZIMUTH-AZIHUTH POSITIONS 
Control Stations: USE MON 1978 and MUSSEL 1932 



Standard Error Used in Computations: 1.3 meters 



Fix 


X 


Y 


Angle of 


Radius of 


No. 


Coordinate 


Coordinate 


Intersection 


90^ circle 


619 


4449.26 


27 1 1. 40 


156. 9 


7.6 


620 


4437. 86 


2807.29 


153.0 


6.5 


621 


4427. 41 


2897. 17 


149.6 


5. 9 


622 


4419.65 


2989.36 


146.2 


5.3 


623 


4410.42 


3072.20 


143. 4 


4.9 


624 


4390.60 


3154.04 


141.3 


4.7 


625 


4376.52 


3234.24 


138.9 


4. 4 


626 


4364.31 


3319. 82 


136.2 


4.2 


627 


4348.71 


341 1. 24 


133.5 


4.0 


628 


4334. 31 


3502. 51 


130.7 


3.3 


631 


4338.97 


3381.55 


135. 1 


4. 1 


632 


4337. 95 


3290. 68 


138.7 


4.4 


633 


4327.70 


3199. 92 


142. 8 


4.8 


634 


4322. 53 


3107.83 


146. 9 


5.4 


635 


4323.18 


3012.49 


151.0 


6. 1 


636 


4324.42 


2916. 00 


155.3 


7.1 


637 


4327.25 


2818.39 


159. 8 


8.7 


638 


4394.77 


2806. 86 


156.0 


7.3 


639 


4386. 58 


2903. 09 


152.0 


6.3 


640 


4377. 29 


2993.63 


148.4 


5.6 


641 


4367. 00 


3090.26 


145. 1 


5. 1 


642 


4355. 55 


3187. 06 


141.9 


4.7 


643 


4345. 97 


3285. 88 


138. 5 


4.4 


644 


4256.90 


3516. 02 


133.6 


4.0 


645 


4260.68 


3416. 65 


137.4 


4.3 


646 


4264. 77 


3321. 01 


141. 1 


4.7 


647 


4283.59 


3208. 82 


144. 8 


5. 1 


648 


4293.24 


3127.36 


147.7 


5.5 


649 


4300.30 


3024.70 


151.8 


6.3 


650 


4345. 53 


3145.58 


144. 1 


5. 1 


651 


4370.38 


3236.09 


139. 1 


4.5 


652 


43S8.77 


3327.77 


134.2 


4.0 


653 


4411.48 


3421. 46 


130. 2 


3.8 


654 


4438.30 


3506. 23 


125.9 


3.5 


655 


4470. 97 


3591.21 


121. 6 


3.3 


656 


4502. 73 


3677. 50 


117.4 


3.2 


657 


4514.38 


3767. 26 


114.0 


3. 1 


658 


4512. 00 


3860.13 


111.1 


3.0 


659 


4520. 10 


3948.53 


107. 9 


2.9 


660 


4494.79 


3049.75 


139.2 


4.5 


661 


4487. 96 


3144. 60 


136. 2 


4.2 


66 2 


4477.62 


3243. 86 


133.2 


4.0 


663 


4462. 90 


3372.74 


129. 4 


3. 7 


664 


4453.74 


3469. 99 


126.5 


3.5 


665 


4440. 56 


3564.89 


123. 8 


3.4 


666 


4465. 02 


3652.95 


119.7 


3.3 


667 


4507. 24 


3743.29 


1 15.0 


3. 1 


668 


4577.55 


3595.27 


116.9 


3.2 


669 


4569.50 


3681.02 


114. 6 


3. 1 


670 


4563.68 


3776.77 


111.9 


3.0 


67 1 


4564.53 


3872. 77 


108.8 


3.0 


672 


4563.90 


3965. 68 


106. 0 


2.9 


673 


4560.41 


4057. 10 


103.2 


2.9 


674 


4553.49 


4147. 06 


100.6 


2.8 


675 


4556.71 


4239.64 


97.7 


2. 8 


677 


4571.60 


4416.22 


92.0 


2.8 


678 


4576.20 


4504.92 


89.3 


2.8 


679 


4582. 23 


4597.03 


86.5 


2.8 


680 


4604.35 


4631. 33 


85. 1 


2. 8 


681 


4613. 39 


4527. 17 


87.9 


2.8 
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Fix 

No. 

682 

683 

684 

685 

686 

687 

688 

689 

690 

69 1 

692 

695 

696 

697 

698 

699 

700 

701 

702 

703 



AZIMUTH-A2IMDTH ACCURACIES (CONTINUED) 

X. Y. Angle of Radius of 



Coordinate 


Coordinate 


Intersection 


90% Cir, 


4611.20 


4421.29 


91.0 


2.8 


4610. 07 


4315. 24 


94.1 


2. 8 


4609.84 


4204. 66 


97. 4 


2.8 


4603. 99 


4098.63 


100.7 


2.3 


4600.52 


3992. 16 


104.0 


2. 9 


4601.27 


3883.40 


107.3 


2.9 


4601. 13 


3780.02 


1 10. 4 


3.0 


4602. 22 


3675.49 


113. 5 


3. 1 


4601.48 


3574. 94 


116.5 


3. 1 


4603. 21 


•3458.09 


120. 0 


3.3 


4524. 84 


3728. 78 


114.8 


3. 1 


4657.06 


3506. 75 


116. 1 


3. 1 


4648. 60 


3602. 24 


1 13. 7 


3. 1 


4630. 15 


3696. 13 


111.8 


3.0 


4629. 80 


3793.39 


108. 9 


3.0 


4623. 03 


3889.32 


106. 3 


2.9 


4622. 83 


3978. 70 


103.7 


2.9 


4617. 60 


407 1. 98 


101 . 1 


2.8 


4623.30 


4163. 73 


98. 2 


2.8 


4618. 23 


4256.25 


95.7 


2. 8 
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CLASSIFIED RANGE-AZIMUTH POSITIONS 
Contrcl Stations: MUSSEL 1932 occupied, initial USE MON 1978 



Standard Errors: Ranee — 3 meters; T-2 — 1.3 meters 



Fix 


X 


Y 


Angle of 


Radius o 


No. 


Coordinate 


Coordinate 


I n^£s^ t i on 


90% Circ 


414 


2898.58 


4611.99 


90.0 


5. 1 


415 


2979.40 


4483.31 


90.0 


5. 1 


416 


3027.69 


4342.47 


90.0 


5. 1 


4 17 


2795.26 


4216.04 


90. 0 


5. 1 


418 


3002.90 


4061. 15 


90.0 


5. 1 


419 


2939. 58 


3938.67 


90.0 


5. 1 


420 


2904.53 


3888.38 


90.0 


5. 1 


421 


3195.47 


3380.99 


90.0 


5. 1 


422 


3128.22 


3487.82 


90. 0 


5. 1 


423 


3048. 17 


3599.79 


90.0 


5. 1 


424 


2994. 68 


3723. 65 


90.0 


5. 1 


425 


2911.61 


3833. 20 


90.0 


5. 1 


426 


2833.62 


3921. 16 


90. 0 


5.1 


427 


2749. 17 


4000.71 


90.0 


5. 1 


428 


2752.46 


4029. 48 


90.0 


5. 1 


429 


2823. 93 


3969.92 


90.0 


5. 1 


430 


2894. 60 


3920. 66 


90. 0 


5. 1 


431 


2907. 95 


3839. 53 


90.0 


5. 1 


432 


2961.44 


3901.64 


90.0 


5. 1 


433 


3005.23 


3980. 1 1 


90.0 


5. 1 


434 


3037.78 


4058.0 1 


90.0 


5. 1 


435 


3065.09 


4134. 51 


90. 0 


5. 1 


436 


3074. 88 


4036. 74 


90.0 


5. 1 


437 


3048. 02 


3953.93 


90.0 


5. 1 


438 


2996.76 


3871. 95 


90. 0 


5. 1 


439 


2943. 10 


3809,76 


90.0 


5. 1 


440 


3128.59 


4049. 87 


90.0 


5. 1 


441 


3096. 20 


3978.78 


90.0 


5. 1 


442 


3058. 84 


3897.24. 


90.0 


5. 1 


443 


3009. 6 1 


3825. 86 


90.0 


5. 1 


444 


2914. 71 


3837.97 


90.0 


5. 1 


445 


3008.50 


3755. 28 


90.0 


5. 1 


446 


3064.52 


3839. 64 


90. 0 


5. 1 


448 


3153.77 


4023.36 


90.0 


5. 1 


449 


3182. 88 


4128.56 


90.0 


5. 1 


450 


3210.38 


4048. 54 


90.0 


5. 1 


451 


3179.45 


3970. 85 


90.0 


5. 1 


452 


3146. 52 


3889.81 


90.0 


5. 1 


453 


3099.79 


3820.34 


90.0 


5. 1 


454 


3047.75 


3751. 18 


90.0 


5. 1 


456 


3037. 99 


3689. 19 


90.0 


5.1 


457 


3121.90 


3785.62 


90.0 


5. 1 


458 


3186.09 


3883. 49 


90.0 


5. 1 


459 


3234.50 


4006. 49 


90.0 


5. 1 


460 


3264. 92 


4128. 69 


90.0 


5. 1 


46 1 


3079.28 


3672. 85 


90.0 


5. 1 


462 


3159.97 


3763. 25 


90.0 


5. 1 


463 


3221.69 


3867.47 


90.0 


5. 1 


464 


3273.65 


3994.32 


90.0 


5. 1 


465 


3324.98 


4005.30 


90.0 


5. 1 


466 


3295.53 


3924.46 


90.0 


5. 1 


467 


3257. 07 


3846. 60 


90.0 


5. 1 


468 


3218. 14 


3770. 87 


90. 0 


5.1 


469 


3172. 68 


3704- 85 


90.0 


5. 1 


470 


3122.00 


3646. 77 


90.0 


5. 1 


47 1 


3097. 64 


3628. 62 


90. 0 


5.1 


472 


3089.39 


3573. 64 


90.0 


5.1 


473 


3118. 69 


3586. 14 


90.0 


5. 1 


474 


3210.37 


3698. 21 


90.0 


5. 1 


475 


3262.22 


3765.90 


90.0 


5. 1 
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Fix 

No. 

U76 

U78 

479 

480 

481 

482 

483 

484 

485 

486 

437 

488 

489 

490 

49 1 

492 

493 

494 

495 

496 

497 



RANGE-AZIfiUTH ACCURACIES (CONTINUED) 



X 


Y 


Angle of 


Radius of 


Coordinate 


Coordinate 


Intersection 


90 % Circle 


3317. 64 


3873. 69 


90.0 


5. 1 


3366. 64 


3900.37 


90.0 


5. 1 


3335.51 


3831. 40 


90.0 


5.1 


3303. 12 


3760.21 


90.0 


5. 1 


3258.28 


3696. 16 


90.0 


5. 1 


3212. 24 


3631.26 


90.0 


5. 1 


3165. 13 


3578.90 


90.0 


5. 1 


3129.45 


3550.77 


90.0 


5. 1 


3142.92 


3499.37 


90.0 


5. 1 


3214. 10 


3575.56 


90.0 


5. 1 


3280. 84 


3665. 79 


90.0 


5. 1 


3344.31 


3753.54 


90.0 


5. 1 


3395.65 


3847. 40 


90.0 


5. 1 


3420.4 1 


3919. 1 1 


90.0 


5. 1 


3466. 94 


3949.20 


90.0 


5. 1 


3439.67 


3856.67 


90.0 


5. 1 


3395.38 


3761. 44 


90.0 


5. 1 


3335.71 


3669.35 


90.0 


5. 1 


3275.38 


3583. 62 


90.0 


5. 1 


3197. 78 


3504. 04 


90.0 


5. 1 


3161.00 


3465. 96 


90.0 


5. 1 
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MPENDIX C 

PBCGBAM FOB 90 PEECENT CONFIDENCE ELLIPSE PAffiAMETEBS 



PECGEAM NAME: ELLIP 

DESCRIPTION: COMPOTES ORIENTATION, 90 % SIGMA-X, 90 % SIGMA-Y, 
FOR ERROR ELLIPSE ABOUT A HYDROGRAPHIC 
POSITION ESTABLISHED BY RANGE-RANGE, AZIMUTH- 
AZIMOTE OR EANGE-AZIMOTH POSITION 

AOTHCE; NICHOLAS E. PEEOGINI 
LT. NOAA 

NAVAL POSTGRADUATE SCHOOL 
DATE: SEPTEMBER , 1984 



IMPLICIT REAL * 4 (A-H,C-Z) 

INITIALIZE VALUES: FOE RANGE-RANGE AND AZIMUTH-AZIMUTH: 

-XL AND YL ARE COORDINATES OF LEFT STATION 
-XR AND YE ARE COORDINATES OF EIGHT STATICN 
-SIGl AND SIGR ARE RESPECTIVE STANDARD ERRORS 
ASSCCIATED WITH EACH LOP 

FOE RANGE-AZIMUTH: 

-XL AND YL ARE COORDINATES OF OCCUPIED STATION 
-SIGL IS SIGMA OF THEODOLITE LOP 
-SIGR IS SIGMA OF RANGE LOP 

******** *********44 **********************■>:********** * 4 :*** 

XL=4914, 75 
Y1=2C09. 86 
SIGL=3. 0 



XR=2474. 75 
YR=4247. 42 
SIGE= 3.0 

ENTEF CORRELATION CCEFFICIENT: USUALLY ZERO FOR R-R, 

R-AZ, AND AZ-AZ 

RO = 0.0 
PI=3. 141593 

ENTER INDICATOR TO TELL WHAT KIND OF DATA IS ENTERING PEOGRA'l 

IND = 1 RANGE-RANGE 

IND = 2 AZIMUTH-AZIMUTH 

IND = 3 RANGE-AZIMUTH 

************************************************** 

IND=1 

************************************************** 



INC IS A TOGGLE WHICH CHECKS FOR BETA GREATER THAN 90 DEG. 
NOTKIKG IN PROGRAM SHOULD BE CHANGED FROM HERE ON 



READ IN DATA FROM 



SENTINEL IS IFIX = 



lATA FILE: IFIX = FIX NUMBER 
PX = X COORDINATE OF HYDRO POSITION 

PY = Y COORDINATE OF HYDRO POSITION 

TD = ANGLE OF INTERSECTION IN DEGREES 
999, TELLS PROGRAM TO STOP READING 



10 CCNTINUE 
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IF (IFIX. EQ. 999) GO TO 900 
IF (ID.LT.90.) LD=TD 
IF flD.LT. 90. ) GO TO 30 

C WORK WITH BETA LESS THAN 90 DEGREES: TOGGLE TURNED CN TC ONE 
DD=180- -TD 
INC=1 

30 CCNTINUE 
C 

C KEEF TANGENT FUNCTICN FROM GOING UNDEFINED IN A RARE CASE 
C OF TEE FIX AND CONTROL STATION HAVING SAME COORDINATES 
IF (EX. EQ- XL) EX=PX+ 0.5 
IF(EY.EQ.YL) EX=PY+ 0.5 

CHANGE DEGREES TO RADIANS 
BETA=.0174532S*DD 

USE LEFT STATION AS BASIS FCR COMPUTATIONS 
ORIENTATION ANGLES WILL BE FIXED WITH RESPECT TO LEFT LOP 




FINE AZIMUTH FROM NORTH BETWEEN HYDRO POSITION AND LEFT 
STATION. AZIMUTH KILL EE DEFINED BETWEEN 0-180 DEGREES 
MEASURED CLOCKWISE FROM NORTH. 

THIS IS THE RANGE-RANGE AZIMUTH DETERMINATION. 

IF (INE.NE. 1) GO TO 40 
IF PY.GE.YL) TEEN 

IF(PX.GE.XL) THEN 

ALPHA = PI-ATAN ( (PY-YL) / (PX-XL) ) 

ELS E 

ALPHA = ATAN ( (PY-YL)/(XL-PX) ) 

END IF 

ELSE 



IF (PX.GE.XL) 
ALPHA = 

ELSE 



THEN 

ATAN ( (YL-PY)/ (PX-XL)) 



ALPHA = 
END IF 
END IF 
GO TO 60 

AZIMUTH FIXING FOE AZI 
40 CONTINUE 



PI-ATAN( (YL-PY)Z(XL-PX)) 



MUTH-AZIMUTH POSITIONS 



IF(PY.GE.YL) TEEN 

IF(PX.GE.XL) THEN 

ALPHA = ATAN ( (PX-XL )/ (PY- YL) ) 

ELSE 

ALPHA = PI-ATAN ( (XL-PX)/ (PY-YL) ) 

END IF 

ELSE 

IF (PX.GE.XL) THEN 

ALPHA = PI-ATAN ( (PX-XL)/ (YL-PY) ) 

ELSE 

ALPHA = ATAN ( (XL-PX) /(YL-PY) ) 

END IF 
END IF 

AZIMUTH EQUALS THETA FOR RANGE AZIMUTH CASE, ASSUMING 
THEODOLITE SIGMA IS LESS THAN RANGE SIGMA 

IF (IND.EQ.3) GO TO 70 



BEGIN COMPUTING THETA, THAT IS THE ANGLE OF ROTATION FROM 
LEFT LOP 
60 CCNTINUE 

B 1=SIGL**2*SIN (2*BETA) +2*RO*SIGL*SIGR*SIN ( BETA) 
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B2=SIGL**2*C0S (2*BETA) ♦2*R0*SIGL*SIGR*C0S (BETA) 4SIGR**2 

IB (ABSiB2) .LI. 0.000 1) E2=.0001 

B3=B1/B2 

C COMfUTE ROTATION ANGLE FROM LEFT LOP 
TH=0.5*ATAN (B3) 

90 CCNTINDE 



DEFINE SEMI-MAJOR AXES ORIENTATION IN TERMS OF 0-180 EEGREES 
ROTATION, CLOCKWISE FROM NORTH 

RANGE-RANGE CASE 

lE(IND.EQ.I) TEEN 

IF(INC.EQ.I) THETA = AIPHA+TH 
IFjlNC.EQ.O) THETA = AIPHA-TH 
END IF 

AZIMUTE-AZIMUTH CASE 

IF(IND.EQ.2) TEEN 

IF(INC.EQ.O) THETA=AIPHA + TH 
IFjlNC.EQ. 1) THETA = AIPHA-TH 
END IF 



FIX RCTATION ANGLE FROM 0-180 DGREES 
70 CCNTIHOE 

C CCNIITICN FOR RANGE-AZIMOT H DATA 
IF{IND.EQ.3) THETA=ALPHA 
IF (THETA. LT.O.) THETA= FI+THETA 
IFjTHETA .GT.PI) THETA= THETA-PI 
C DEG IS THE SEMI-MAJOR ELLIPSE AXIS ORIENTATION IN DEGREES 
DEG=57. 295779*THETA 

C COMEUTE 90% SIGMAX AND SIGMAY OF ERROR ELLIPSE 

CALL PROB (SIGI,SIGR,RO,TD,SGX90,SGY90,RADIUS,ELAR,CIRAR) 



100 



900 



WEITE(7.ldO) iriX,^X.I>Y,fD.SGX§0,SGY§0,DEG 
FCRMATnX,I3,3X,F7.2,3X,F7.2,3X,F5. 1,3X,F4. 1, 



3X,F5. 1) 
GO TO ID 
CONTINUE 
STOP 
END 



3X,F4.1 
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APPENDIX D 

ACCORACY CLASSIFICAIION: 90 PERCENT CONFIDENCE ELLIPSES 



CLASSIEIED RANGE-RANGE POSITIONS 
Control Stations: BEACH LAB 1982 and MUSSEL 1932 

Standard Error Used in Computations: 3 meters 



Fix 


X 


Y 


Beta 


90% 


90% 


Ori€nta- 


No. 


Coord. 


Coord. 




Siqma X 


Siqma Y 


tion 


1 


2668.05 


4942.07 


127.0 


10.2 


5. 1 


79.0 


2 


2852.01 


5076.69 


121.6 


9.3 


5.2 


85.3 


3 


3041.58 


5194.85 


1 18.6 


8.9 


5.3 


90.2 


a 


3040.27 


5148.74 


1 17.0 


8.7 


5.3 


90.6 


5 


2838.26 


5021.35 


120.3 


9. 1 


5.2 


85.3 


6 


2640.06 


4877. 13 


126.9 


10.2 


5. 1 


78.1 


7 


2553.46 


4752.36 


130.4 


10.9 


5.0 


74.1 


8 


2724.10 


4885.47 


121.4 


9.3 


5.2 


82.0 


9 


2889.06 


5013.63 


1 17. 6 


8.8 


5.3 


87.2 


10 


3C75.61 


5124.37 


1 15.0 


8.5 


5.4 


91.9 


1 1 


3172.63 


5 136.23 


1 12.7 


8.2 


5.5 


94.5 


12 


2958. 13 


5006.65 


1 14.4 


8.4 


5.4 


89.7 


13 


2771.33 


4876.02 


1 18.0 


8.8 


5.3 


84.2 


14 


2581.02 


4729.89 


127.0 


10.2 


5. 1 


75.9 


15 


2584.41 


4665.31 


124.0 


9.7 


5.2 


76.7 


16 


2740.05 


4805.63 


1 16. 7 


8.7 


5.3 


83.8 


17 


2913.27 


4928.51 


1 12.8 


8.2 


5.5 


89.2 


18 


3097.44 


5047.77 


111.2 


8. 1 


5.5 


93.5 


19 


3193.94 


5 103.46 


1 10.9 


8.0 


5.5 


95.5 


20 


2627.47 


4659.05 


• 118.8 


8.9 


5.3 


79.8 


2 1 


2762.15 


4769.60 


1 13.2 


8.3 


5.5 


85.4 


22 


2904.91 


4875.24 


1 10.5 


8.0 


5.5 


89.7 


23 


3056.38 


4975.73 


109.3 


7.9 


5.6 


93.3 


24 


3207.36 


5064.00 


108.9 


7.8 


5.6 


96.3 


25 


3373.36 


510 1.96 


107.1 


7.7 


5.7 


100.0 


26 


3190.64 


5007.38 


1 06.8 


7.6 


5.7 


96.7 


27 


3015.45 


4900.75 


107.1 


7.7 


5.7 


93.2 


28 


2839.38 


4778.31 


1 08.7 


7.8 


5.6 


88.8 


29 


2679.99 


4651.31 


1 12.8 


8.2 


5.5 


83.4 


30 


2727.88 


4632.77 


106.9 


7.6 


5.7 


86.7 


31 


2899.83 


4776.97 


105.2 


7.5 


5.7 


91.3 


32 


3092.77 


490 1.75 


104.4 


7.4 


5.8 


95.6 


33 


3295.42 


5019.92 


105. 0 


7.5 


5.7 


99.2 


34 


3502. 07 


5123.36 


106.0 


7.6 


5.7 


102.6 


35 


3697.02 


515 1.09 


1 05.3 


7.5 


5.7 


106.2 


36 


3474.77 


5065.57 


104.1 


7.4 


5.8 


102.7 


37 


3257.48 


4953.37 


102.7 


7.3 


5. 8 


99.3 


38 


3043.03 


482 1.88 


101.7 


7. 2 


5.9 


95.5 


39 


2845.07 


4680.80 


101.7 


7.2 


5.9 


9 1.4 


40 


2746.96 


4608.32 


103.1 


7.3 


5.8 


88.6 
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CLASSIFIEE AZiaUTH-AZIMDTH POSITIONS 
Control Stations: USE MON 1978 and MUSSEL 1932 

Standard Error Used in Computations: 1.3 meters 



Fix 


X 


X 


Beta 


90 % 


90 % 


Orienta- 


No. 


Ccord. 


Coord. 




Sijgma X 


Sigma Y 


tion 


619 


4449.26 


271 1.40 


156.9 


9.8 


2.0 


139. 4 


620 


4437.86 


2807.29 


153.0 


8.5 


2.0 


139.8 


62 1 


4427.41 


2897. 17 


1 49.6 


7.5 


2.0 


139. 8 


622 


4419.65 


2989.36 


146.2 


6.8 


2. 1 


139.8 


623 


441 C. 42 


3072.20 


143.4 


6.3 


2. 1 


139. 6 


624 


4390.60 


3 154.04 


141.3 


5.9 


2. 1 


139. 1 


625 


4376. 52 


3234.24 


138.9 


5.6 


2. 1 


138. 6 


626 


4364.31 


331 9.82 


136.2 


5.3 


2.1 


138.0 


627 


4348.71 


341 1.24 


133.5 


5.0 


2. 1 


137. 3 


628 


4334.31 


3502.51 


130. 7 


4.7 


2.2 


136.5 


63 1 


4338.97 


3381.55 


135. 1 


5.2 


2. 1 


137.4 


63 2 


4337.95 


3290.68 


138.7 


5.6 


2. 1 


137. 8 


633 


4327.70 


3199.92 


1 42.8 


6.2 


2.1 


138. 1 


634 


4322.53 


3107.83 


146.9 


6.9 


2. 1 


138.2 


635 


4323.18 


3012.49 


151.0 


7.9 


2.0 


138. 3 


636 


4324.42 


2916.00 


155.3 


9.2 


2.0 


138. 1 


637 


4327.25 


2818.39 


159.8 


11.3 


2.0 


137.7 


638 


4394.77 


2806.86 


156.0 


9.5 


2.0 


138.9 


639 


4386.58 


2903.09 


1 52.0 


8.2 


2.0 


139. 1 


640 


4377.29 


2998.63 


148.4 


7.2 


2. 1 


139. 1 


641 


4367.00 


3090.26 


1 45. 1 


6.6 


2.1 


138. 9 


642 


4355.55 


3187.06 


141.9 


6.0 


2. 1 


138.5 


643 


4345.97 


3285.88 


138.5 


5.6 


2. 1 


138.0 


644 


4256.90 


3516.02 


1 33.6 


5.0 


2. 1 


135.5 


645 


4260.68 


3416.65 


137.4 


5.4 


2.1 


136. 2 


646 


4264.77 


332 1.01 


141. 1 


5.9 


2.1 


136.8 


647 


4283. 59 


3208.82 


144.8 


6.5 


2. 1 


1 37. 5 


648 


4293.24 


3127.36 


147.7 


7.1 


2. 1 


137.8 


649 


4300.30 


3024.70 


151.8 


8.1 


2.0 


137.9 


650 


4345.53 


3145.58 


144. 1 


6. 4 


2. 1 


138.5 


651 


4370.38 


3236.09 


139.1 


5.7 


2. 1 


1 38. 5 


652 


4398.77 


3327.77 


134.2 


5. 1 


2.1 


138. 4 


653 


441 1.48 


342 1.46 


130.2 


4.7 


2.2 


138.0 


654 


4438.30 


3506.23 


125.9 


4.3 


2.2 


137.7 


655 


4470. 97 


359 1.21 


121.6 


4.0 


2.3 


137.4 


656 


4502.73 


3677.50 


1 17.4 


3.8 


2.3 


137.0 


657 


4514.38 


3767.26 


114.0 


3.6 


2.4 


136. 2 


658 


4512.00 


3860.13 


111.1 


3. 5 


2.4 


135.2 


659 


4520. 10 


3948.53 


107.9 


3.4 


2.4 


134.3 


660 


44S4.79 


3049.75 


139.2 


5.7 


2.1 


141. C 
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CLASSIFIED RANGE-AZIMUTH POSITIONS 



Control Stations 


: USE MON 


1 978, 


MUSSEL 


1932 




Standard Errors: 


RANGE — 3 


me ters ; 


T-2 — 1 


.3 




Fix 


X 


Y 


Seta 


9 0% 


90% 


Or lent a- 


No. 


Coord. 


Coord. 




Sigma X 


Sigma Y 


tion 


414 


2898.58 


461 1.99 


90.0 


6. 4 


2. 8 


49.3 


415 


2979.40 


4483.31 


90.0 


6.4 


2.8 


64.9 


416 


3027.69 


4342.47 


90.0 


6.4 


2.8 


80. 2 


417 


2795.26 


4216.04 


90.0 


6.4 


2.8 


95. 6 


418 


3CC2.90 


406 1.15 


90.0 


6. 4 


2.8 


1C9.4 


4 19 


2939.58 


3938.67 


90.0 


6.4 


2.8 


123.6 


420 


2904.53 


3888.38 


90.0 


6.4 


2.8 


129.9 


421 


3195.47 


3380.99 


90.0 


6.4 


2.8 


140. 2 


422 


3128.22 


3487.82 


90.0 


6. 4 


2.3 


139.3 


423 


3048. 17 


3599.79 


90.0 


6.4 


2.8 


138. 5 


424 


2994.68 


3723.65 


90.0 


6. 4 


2.8 


135. 2 


425 


2911.61 


3833.20 


90.0 


6.4 


2.8 


133. 5 


426 


2833.62 


3 921.16 


90.0 


6. 4 


2.8 


132. 3 


427 


2749. 17 


400C.71 


90.0 


6.4 


2.8 


132.0 


428 


2752.46 


4029.48 


90.0 


6. 4 


2. 8 


128. 1 


429 


2823.93 


3969.92 


90.0 


6.4 


2.8 


128. 5 


430 


2894.60 


3920.66 


90.0 


6. 4 


2.8 


127.9 


431 


2907.95 


3839.53 


90.0 


6. 4 


2.8 


133. 3 


43 2 


296 1. 44 


390 1.64 


90.0 


6. 4 


2. 8 


125.4 


43 3 


3C05.23 


3980.11 


90.0 


6.4 


2.8 


116.7 


434 


3C37.78 


4058.01 


90.0 


6. 4 


2.8 


103.6 


435 


3C65.09 


4134.51 


90.0 


6.4 


2.8 


100. 8 


43 6 


3C74.88 


4036.74 


90.0 


6.4 


2.8 


109.3 


437 


3C48.02 


3953.93 


90.0 


6.4 


2.8 


1 17. 1 


438 


2996.76 


3871.95 


90.0 


6. 4 


2.8 


125. 7 


439 


2943. 10 


380 9.76 


90.0 


6. 4 


2.8 


133. 1 


440 


3128.59 


4049.87 


90.0 


6. 4 


2.8 


106. 8 


44 1 


3096.20 


3978.78 


90.0 


6.4 


2.8 


113. 4 


442 


3C58.84 


3897.24 


90.0 


6. 4 


2.8 


120.9 


44 3 


3009.61 


3825.86 


90.0 


6.4 


2.8 


128.2 


444 


2914.71 


3 837. 97 


90.0 


6. 4 


2.8 


132.9 


445 


30C8.50 


3755.28 


90.0 


6.4 


2.8 


132. 7 


446 


3C64.52 


3839.64 


90.0 


6. 4 


2.8 


124.7 


448 


3 153.77 


4023.36 


90.0 


6.4 


2.8 


1 08.3 


449 


3182.88 


4 128.56 


90.0 


6.4 


2.8 


99. 5 


450 


3210.38 


4048.54 


90.0 


6. 4 


2.8 


105. 1 


451 


3 179. 45 


3970.85 


90.0 


6. 4 


2.8 


111.4 


452 


3 146.52 


3889.81 


90.0 


6.4 


2.8 


1 18.0 


453 


3C99.79 


3820.34 


90.0 


6. 4 


2.8 


124.3 


454 


3047.75 


3751. 18 


90.0 


6. 4 


2.8 


130.9 
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